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SECTION  1.0 


INTRODUCTION 


The  conventional  methods  of  performing  longitudinal  structure  designs 
of  ships  make  use  of  accumulated  experience  from  previously  built  ships  of 
similar  size  and  function.  The  accumulated  experience  is  mostly  expressed  in 
the  form  of  semi -empi ri cal  formulas  contained  in  classification  society  rules 
and  design  specifications.  The  designs  resulting  from  this  approach  are  uncertain 
as  to  the  degree  of  structural  adequacy  they  afford  even  though  the  ship  designs 
based  on  these  approaches  have  given  acceptable  service.  The  uncertainty  stems 
from  the  assumptions  made  regarding  parameters  affecting  the  environment  and 
the  strength  of  the  ship.  Many  years  of  design  experience  have  shown  that  by 
using  appropriate  empirical  margins  for  strength  over  expected  load,  the  unknowns 
can  be  accounted  tor  and  ships  with  acceptable  risk  or  probability  of  failure 
levels  designed. 

With  the  advent  of  new  ship  types,  and  the  resultant  lack  of  "accumulated 
experience"  on  vessels  of  similar  size  and  function,  it  has  become  a  professional 
responsibility  to  look  Into  a  more  scientific,  or  rational,  approach 
to  longitudinal  strength  design  of  ship  hulls.  In  this  context,  various  invest¬ 
igators  in  the  ship  research  community  have  adopted  probabilistic  structural 
analysis  procedures  from  mechanical  and  civil  engineering.  In  the  "probabilistic 
approach",  since  the  quantitative  values  of  many  of  the  factors  affecting  the 
strength  of  the  structure  and  the  magnitude  of  the  load  are  statistically 
determined,  the  resulting  measure  of  the  adequacy  of  the  design  is  also 
statistical  in  nature. 

In  the  study  presented  in  this  report,  various  facets  of  probabilistic 
structural  design  were  investigated  with  emphasis  on  applicability  to  ships. 

Section  2.0  g'*ves  a  statement  concerning  the  detailed  objectives  of  the 
study.  In  Section  3.0,  probabilistic  structural  analysis  is  reviewed  from  c  general 
standpoint  and  Its  applicability  to  ships  is  noted.  Section  *1.0  discusses  the 
possible  structural  modes  of  failure  of  a  ship  that  pertain  to  longitudinal 
strength.  The  present  situation  with  information  on  ship  loads  as  they  relate 
to  structural  design  is  discussed  in  Section  5.0,  and  the  probabilistic 
structural  analysis  procedures  that  show  promise  for  ship  applications  are 
presented  in  Section  6.0.  In  Section  7-0,  the  investigations,  analyses  and 
collected  information  performed  and  obtained  as  part  of  this  study  :n  the  area 
cf  the  uncertainties  of  hull  strength  with  respect  to  the  statistical  description 
of  the  strength  are  presented.  Section  8.0  gives  sample  calculations  for 
different  ships  using  a  probabilistic  structural  analysis  procedure  embodied 
In  a  computer  program  included  in  the  Appendix.  Sections  9.0  and  !u  0 
give  the  conclusions  and  recommendations  respectively  arrived  at  as  a  result  of 
these  studies. 

The  references  cited  in  the  report  are  listed  in  Section  11.0 
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SECTION  2.0 


STATEMENT  OF  OBJECTIVES 


The  objectives  pf  this  study  were  modified  by  the  Ship  Structure  Com¬ 
mittee  during  the  course  of  the  project  to  be  commensurate  with  what  was  found 
to  be  available  and  possible  within  the  rather  small  funding  allocated. 

The  final  objectives  can  be  stated  as  follows: 

o  Survey  the  existing  literature  on  reliability  analysis  and  proba¬ 
bilistic  design  methods  in  structures.  Comment  on  the  applicabi¬ 
lity  to  ships. 

o  Develop  a  method,  or  use  an  existing  method,  for  the  formulation 
of  strength  in  terms  of  the  means  and  variances  of  its  uncertain¬ 
ties.  Although  a  mathematical  distribution  of  strength  is  not  re¬ 
quired,  observations  are  to  be  made  with  respect  to  the  impact  of 
using  only  means  and  variances. 

o  Relate  the  existing  bending  moment  distributions  calculated  from 

existing  data  to  the  developed  strength  distributions  using  an  existing 
method  for  structural  reliability  analysis.  Use  available  statis¬ 
tical  strength  parameter  means  and  variances  and  make  assumptions 
for  any  strength  or  load  parameters  for  which  no  statistical  data 
are  avai table. 

o  Develop  a  FORTRAN  IV  computer  program  to  perform  the  above  proce¬ 
dure  with  the  objective  of  determining  the  safety  level  of  a  given 
ship  subjected  to  a  given  load. 

o  Apply  the  developed  computerized  procedure  to  actual  ships. 

o  On  the  basis  of  obtained  results,  suggest  further  research  to 

develop  suitable  longitudinal  strength  criteria  for  future  designs. 
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SECTION  3-0 


PROBABILISTIC  APPROACH  TO  STRUCTURAL  DESIGN 


3- 1  Genera  1 


The  objectives  of  this  study  include  the  analysis  of  uncertainties 
associated  with  ship  hull  strength  and  the  development  of  expressions  for 
structural  reliability.  Such  analyses  require  the  adoption  of  a  probabilistic 
structural  design  approach  since  a  purely  deterministic  approach  cannot  yield 
the  desired  information. 

In  the  deterministic  design  of  structures,  the  strength  of  the  structure 
is  always  increased  above  that  which  would  just  survive  the  greatest  expected 
load  by  an  empirical  margin.  The  ratio  of  the  latter  to  the  former  strength 
is  usually  termed  the  factor  of  safety.  It  accounts  for  all  the  unknowns  in 
the  load  and  strength  and  yields  a  structure  that  should  have  an  acceptable 
performance  based  on  past  experiences. 

The  fundamental  aims  of  a  probabilistic  approach  are  to  more  clearly 
and  rationally  define  the  necessary  margin,  or  factor  of  safety,  and  obtain  a 
quantitative  measure  of  performance  through  a  rational  rather  than  empirical 
analysis.  The  measure  of  performance  Is  usually  called  the  probability  of 
failure  or  reliability.  With  such  aims,  it  is  not  necessary  that  a  probabi¬ 
listic  analysis  be  exhaustive  In  that  rational ization  of  even  only  one  of  the 
unknowns  in  the  factor  of  safety  will  put  i t  on  a  sounder  footing.  In  this 
vein,  the  ultimate  result  of  Improved  probabilistic  analysis  procedures,  as 
far  as  designers  are  concerned,  wi ) 1  probably  be  rat!onal  factors  of  safety 
based  on  desired  quantitative  levels  of  performance.  The  probabilistic 
analysis  itself  need  not  be  executed  by  the  designers,  although  this  could 
be  possible. 

A  complete  probabilistic  structural  analysis  would  proceed  in  the 
fol lowing  manner  [9]*: 

0  Conduct  an  analysis  of  failure  modes,  effects,  and  criticality. 

-  Identify  all  significant  failure  modes  of  the  structure. 

-  List  the  cause  of  these  failure  modes. 

-  Identify  all  parameters  contributing  to  these  causes. 

-  Determine  the  criticality  of  all  siginficant  failure  modes 
to  the  success  of  structures. 

-  List  the  most  critical  failure  nudes  in  order  of  priority. 

0  Formulate  the  relationship  between  the  critical  parameters  and 

the  failure-governing  criteria  involved. 

*  Determine  the  failure-governing  load  function. 

*  Determine  the  failure-governing  load  distribution. 

0  Determine  the  failure-governing  strength  function. 

0  Determine  the  failure-governing  strength  distribution. 

*  Calculate  the  probability  of  failure  or  reliability  associated  wi th 


*  Numbers  in  brackets  Indicate  similarly  numbered  references  in  Section  11.0. 
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the  failure-governing  load  and  strength  distribution  for  each 
critical  failure  mode. 

°  An  upper  bound  cf  the  total  probability  of  failure  or  a  lower 
bound  of  the  reliability  will  be  the  sum  of  the  individual 
probabilities  of  each  of  the  critical  failure  modes  under  the 
assumption  that  these  modes  are  mutually  exclusive  events. 

Because  of  the  difficulty  associated  with  the  determination  of 
the  failure-governing  load  and  strength  functions  and  d i s tr i but  ions, a  number 
of  probabi 1 ist ic  approaches  or  methods  have  evolved.  They  differ  fundamentally 
in  tne  two  primary  aims  of  any  probabilistic  analysis  as  mentioned  above: 

°  Quantitative  measure  of  performance 

°  Rational  quantification  of  load  and  strength 

Actually,  not  all  the  approaches  are  necessarily  probabilistic  in 
the  mathematical  sense  in  that  for  some,  probability  densities  and  distributions 
are  not  needed,  and  the  output  is  not  a  probability. 

These  methods  may  be  grouped  as  follows: 

0  Classical  probabi 1 istlc  approach 

°  Safety  index  approach 

0  Strength  reduction  and  load  magnification  factors  approach 

The  presentation  in  this  section  is  divided  into  three  groups.  The 
first  group  discusses  the  general  approach  used  in  obtaining  the  quantitative 
measure  of  performance  of  a  structure  given  the  load  and  strength  statistics. 

The  next  groups  each  deal  with  details  of  the  strength  and  load  formulations 
respectively,  in  a  general  sense.  More  specific  mention  of  these  considerations, 
as  applicable  to  ships,  is  given  in  Sections  5-0  thru  7.0,  respectively  for 
leadings,  longitudinal  strength,  and  for  uncertainties  in  the  strength  of 
the  ship's  hui 5 . 

The  literature  contains  abundant  sources  of  probabilistic  structural 
analyses.  Most  of  the  work  has  been  done  in  the  areas  of  civil  and  mechanical 
engineering  but  has  more  recently  spread  to  naval  architecture. 

Probabilistic  design  concepts  for  structures  were  first  proposed  In 
the  U.S.  in  1947  t 13  .  Since  then,  several  investigators  have  presented 
further  considerations  for  applications  in  civil  engineering.  References  [2] 
thru  [6],  mechanical  engineering,  references  [7]  thru  [9],  and  more  recently 
in  naval  architecture,  reference  1 10]. 

Within  the  framework  of  the  present  study,  a  brief  review  of  the 
numerous  methods  as  cited  was  performed  to  identify  the  ones  which  would  seem 
appropriate  for  future  consideration  in  probabilistic  structural  analyses 
of  ships  from  the  standpoint  of  design. 

3.2  Probabi I istic  Methods 

3.2.1  Quantitative  Measure  of  Performance 

As  previously  mentioned,  the  existing  probabilistic  structural 
analysis  methods  differ  in  the  output  measure  of  performance  of  the  structure 


i 


being  considered. 


Those  methods  that  are  more  probabilistic  in  the  mathematical 
s^nse,  generally,  are  of  the  classical  type.  Their  measure  of  performance  is 
in  terms  of  a  probability  defining  failure  or  reliability. 

The  other  methods  have  evolved  primarily  due  to  the  difficulties 
associated  with  executing  a  fully  probabilistic  procedure.  Their  measure  of 
performance  is  not  a  probability  at  alt,  instead, it  is  a  number  Indicating 
ei »er  a  margin  of  safety  or  reduction  and  magnification  factors  for  strength 
and  load,  respectively.  These  numbers  do  not  have  a  physical  significance 
like  probability  of  failure  or  reliability,  but  they  can  be  compared  to  each 
other  for  previous  successful  and  unsuccessful  designs  to  obtain  limiting  values. 

3.2.2  Classical  Approach 

The  one  common  point  in  all  probabilistic  structural  analysis  pro¬ 
cedures  is  the  definition  of  the  probability  of  failure  and  reliability.  If 
the  fa i lure -governing  load  is  2  and  the  failure-governing  strength  S,  then  the 
probability  of  failure,  Pf,  is  given  by  all  probabilities  that  the  failure- 
governing  load  exceeds  the  failure-governing  strength: 

Pf  -  P  (2>S)  ^ 

The  probability  of  failure  is  also  called  the  unreliability,  while 
the  reliability,  R,  becomes: 

R  ■  1-Pf  -  P  (S>2)  (2) 

Equation  (I)  is  presented  in  much  of  the  literature,  for  example 
in  [10],  as  directly  applicable  to  ships  in  the  following  manner; 

Pf  -  P  [$<2]  -  P  [J<!]  -  P  [Q< 1 ]  (3) 

-  P  [(S-2)<0]  -  P  [M<0] 

The  terms  nQV  and  "M"  of  Equation  (3)  are  functions  of  two  random  variables: 

the  strength,  S,  and  the  load,  2,  and  themselves  random  variables  whose 

probability  must  be  determined  by  joint  probability  density  and  distribution 

functions.  However,  there  seems  to  be  a  universal  agreement  to  consider 

the  load  and  strength  statistically  independent  so  tnat  the  statistics  of  M  and  Q  can 

be  directly  determined  from  those  of  S  and  2.  This  assumption  appears  to  be 
reasonable  for  most  strength  considerations  as  long  as  the  effects  on  the 
structure  of  being  in  an  aqueous  environment  with  waves  for  a  long  period  of 
time  are  accounted  for  in  the  strength.  If  4>  (z)  and  ClO  are  the  probability 
density  and  distribution  functions  of  the  load,  respectively,  and  fc  (s)  and  F  (s) 
those  of  strength,  then  it  can  be  shown  that  the  density  and  distribution  s 
functions  of  Q  are,  [10]: 

fQ  (q)  “7*  (2)  f«  (qz)  z  dz  (*») 

V.  0 

Fq  (q)  “  0/*(z)  Fs  (qz)  dz  (5) 
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(6) 

(7) 


and  the  probability  of  failure  becomes: 

Pf  •  ?♦(*)  Fs  (z)  dz 

-  I  -J*U)  fs  (z)  d z 

Equations  (6)  and  (7)  are  rather  simple  and  could  easily  be 
evaluated  provided  the  density  and  distribution  functions  of  load  and 
strength  are  known.  This  Is  where  the  crux  of  the  matter  lies  and  will 
ba  discussed  later  in  Sections  3.3  and  3-i>-  The  methods  that  make  use 
of  Equations  (6)  and  (7)  vary  significantly  in  complexity  and  effort 
required  for  execution. 

Equation  (7)  can  be  evaluated  for  each  mode  of  failure  and, 
as  noted  previously,  the  sum  of  all  probabilities  of  failure  for  all  modes 
will  give  an  upper  bound.  To  do  better  would  require  the  joint  probability 
density  function  of  strength  in  the  various  failure  modes  which  would  be 
at  best  very  difficult  to  obtain.  A  lower  bound  on  the  probability  of 
failure  can  be  determined  by  assuming  that  the  modes  of  failure  are  perfectly 
correlated. 

3-2.3  Safety  Index  Approach 

The  difficulty  in  obtaining  load  and  strength  density  and  dis~ 
tribution  functions  has  led  investigators  to  develop  approaches  which  mini¬ 
mize  the  effort  required.  For  instance,  in  the  area  of  ships,  [13]  contains 
an  approximate  semi -probabi I i st i c  design  method  which  was  motivated,  among 
other  things,  by  the  lack  of  data  on  loads  and  strength  and  by  the  contro¬ 
versial  status  of  forms  of  load  and  strength  distributions.  The  method 
requires  that  only  the  means  and  variances  of  the  load  and  strength  be  known. 

This  "approximate"  approach  considers  the  margin  of  safety  M 
of  Equation  (3)  as  a  random  variable  with  mean  nyj  and  variance  c^. 

Pf  -  P  lM<0]  -  P  [^^?<^EL]  -  P  [g*-Y]  -  Fc  (- Y  )  (8) 

By  using  the  error  distribution  of  M,  [16],  discussed  in  more  detail  in 
Section  3>3>  the  mean  and  the  variance  of  H  can  be  written: 

"V,  *  mj  -  m2  (9) 

o„2  *  o*  +  o*  (10) 

where: 

m$,  Oj-  mean  and  variance,  respectively  of  strength. 

mE *  aZ  m  mean  and  variance  of  total  load. 

The  following  results  are  obtained  by  algebraic  processes: 


"H  _  "S  ~«"Z  .  _ e-1  (11) 

"H  Y'S2+C2?  "  /o24+v22 
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02) 


)|/2 


l-rvs 

mj  »e  mg 

(13) 

m„ 

SM  -  -JL 

0*0 

°N 

safety  index  ■  m^/i^ 

central  safety  factor  -  m^/m^ 

o« 

Vc  *  coefficient  of  variation  (CO V)  of  strength  =■  — — 

*  m5 

SM  -  required  section  modulus  of  the  ship  hull 

V,  ■  COV  of  load  »  !QL 
r  mg 

SM  -  required  section  moddlus  of  ship  hull 
cr^  “  average  of  failure  stress  of  hull  material 

-  mean  of  the  margin  of  safety 
a*  -  variance  of  the  margin  of  safety 

From  Equation  (8)yit  can  be  seen  that  each  value  of  the  safety 
index  Y  is  associated  with  some  probability  of  failure.  However,  Equation  (8) 
cannot  be  evaluated  since  the  distribution  function  is  not  known.  If 
enough  information  were  available  to  determine  Fp,  then  Equations  (61  and  (7) 
of  the  classical  approach  could  be  used  directly.  From  Equations  (II) 
through  (lM,  It  can  be  seen  that  the  inputs  needed  to  obtain  a  hull  design 
strength  are  the  strength  and  load  COV's,  mean  of  the  bending  moment,  and 
the  safety  index  Y*  The  amount  of  computation  j5  insignificant. 

The  safety  Index  Y  Is  a  sinqle  number  that  must  be  obtained 
On  the  basis  of  many  technical  factors.  It  has  previously  been 
proposed  ( 1 3 1  to  determine  this  value  from  existing  designs  to  take  into 
account  the  vast  accumulated  experience.  In  addition,  if  the  probability 
of  failure  associated  with  past  designs  is  socially  acceptable,  then  this 
aspect  is  also  considered. 

3 . 2 . h  Strength  Reduction  and  Load  Magnification  Factors 

This  method,  discussed  in  [5,62,63],  is  similar  to  the  approximate 
method  described  above  In  that  only  means  and  variances  of  the  load  and 
strength  are  used  to  obtain  relative  and  serni -probabi 1 istic  measures  of  the 
structure's  performance.  In  this  case,  the  measures  of  performance  are  the 
strength  reduction  and  load  magnification  factors. 
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The  strength  reduction  factor, 
fg,  can  be  defined  as  follows: 


and  load  magnification  factor. 


f 


S 


minimum  strength  _  mS~^S~  °S 
average  strength  m 


I  -  KSVS;  fs<  I 


(15) 


_  maximum  load 
r  average  load 


mz+Kz  z 


I  + 


K  V  • 

■zV 


(16) 


where:  =  Factors  giving  the  number  of  standard 

deviations  between  the  average  and  the 
minimum  strengths  and  the  maximum  loads, 
respectively. 

For  a  safe  design,  the  minimum  strength  must  exceed  or  equal 
the  maximum  load: 


msfs  >  m2  fz  (17) 

The  values  of  acceptable  strength  reductions  factors  and  load 
magnification  factors  could  be  obtained  from  past  designs  in  a  similar 
fashion  to  the  safety  index  of  the  previous  section. 

In  [5],  this  approach  has  been  extended  to  fatigue  for  both  the 
constant  range  and  the  random  loads. 

Similarly  to  the  safety  index  approach,  the  analyses  required 
to  execute  this  method  are  quite  limited  in  extent  and  complexity. 

3 . 3  Strength  Statistics 

3.3.1  General 

It  must  be  first  stated  that  the  strength  of  the  hull  girder 
may  or  may  not  vary  with  time  depending  on  the  failure  mode  being  con¬ 
sidered.  Time  invariant  strengths  will  include  yielding  and  buckling. 

Time  variant  strengths  will  include  fracture,  fatigue,  and  reduced  strengths 
due  to  corrosion.  For  ships,  time  variant  strengths  will  also  normally 
include  random  loadings  of  low  or  high  cycles,  and  possibly  thermal  loadings. 
This  scenario  should  cover  the  most  significant  modes  of  hull  girder  failure 
which  need  to  be  addressed. 

3.3-2  Strength  Equation 

The  strength  of  a  structure  is  principally  described  in  two 
different  ways  in  the  numerous  probabilistic  structural  design  methods  to 
be  found  in  rhe  literature. 


S  -  f  <V  — -  *„) 

(18) 

S-  kjkjkj— -kns' 

(19) 
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where: 


Constituent  parts  of  the  strength  which 
are  assumed  to  be  random  variables 


S  *  Nominal  strength  determined  under  idealized 
and  standard  test  conditions 

Kj - 1^  “  Strength  factors  to  convert  the  nominal 

strength  to  actual  strength.  (These  factors 
are  assumed  to  be  random  variables). 

The  K  factors  account  for  physical  variables  such  as  size,  forming  and  manu¬ 
facturing  processes,  surface  finish,  load,  heat  treatment,  direct  surface 
environment,  temperature,  time,  corrosion,  etc. 

The  approach  given  by  Equation  ( 1 8)  has  been  used  in  ships, 
but  the  actual  examples  developed  have  been  such  that  only  the  explicit 
functional  strength  constituents,  S,  have  been  considered  as  random 
variables  or  uncertainties  In  the  strength.  As  the  probabilistic  analyses 
become  more  comprehensive  and  more  uncertainties  become  identified,  some 
of  these  may  not  appear  as  constituents  In  the  strength  equation,  and 
the  approach  depicted  in  Equation  (19)  may  have  to  be  adopted  in  addition 
to  that  in  Equation  ( 1 8) . 

3*3.3  Strength  Distributions 

Equations  (18)  and  (19)  give  general  expressions  for  the 
strength,  but  since  the  strength  is  statistical  in  nature,  the  probability 
density  and  distribution  function  must  be  specified  to  completely  characterize 
it  and  allow  the  probability  of  failure  to  be  evaluated  by  Equations  (G)  and 
(7). 


The  probabilistic  structural  analysis  approaches  found  in 
the  literature  assume  that  the  strength  distribution  can  be  determined  in 
one  of  the  following  ways: 

*  Actual  component  strength  distribution  determined  by 
actual  testing  under  the  exact  geometry,  application, 
and  operational  environment  in  which  the  component 
shall  function. 

*  Component  strength  distribution  synthesized  from  the 
known  distributions  of  the  constituent  parts  and 
strength  factors  as  given  in  Equations  (18)  and  (19). 

*  An  assumption  made  as  to  what  type  of  distribution  the 
strength  will  follow,  i.e.  normal,  lognormal,  Weibull,  etc. 

*  An  assumption  made  that  all  that  can  be  determined  of 
the  strength  is  its  COV. 

The  first  of  the  above  approaches  is  used  extensively  in 
machine  design  and  some  of  the  test  equipment  required  is  described  in  [7). 

This  approach  would  hardly  seem  realistic  for  ships  because  of  the  large 
size  of  the  structure,  the  implication  of  using  the  whole  ship  as  a  dis- 


cardable  test  component,  and  the  large  data  sample  required  for  conclusive 
results.  Whether  or  not  components  of  the  ship  structure  could  be  tested 
and  results  extrapolated  to  the  whole  ship  appears  questionable.  In  the 
case  of  welded  ship  grillages  under  compressive  load  [64]: 


"Further  experimental  evaluation  of  grillage  strength  also 
has  a  key  part  to  play  but  cannot  be  expected  to  provide  direct  statistical 
descriptions  of  grillage  strength;  large-scale  tests  of  the  type  described 
in  the  present  paper  are  too  expensive  to  carry  out  in  sufficient  numbers 
and  small-scale  tests  are  statistically  unrepresentative  for  the  reasons 
mentioned  above.  It  is  suggested  that  the  main  role  of  further 
grillage  tests  should, therefore, be  to  guide  the  development  of  improved 
analysis  methods  and  to  check  the  accuracy  of  such  methods  and  design 
data  with  provision  of  empirical  corrections  where  necessary." 

The  second  approach  requires  that  the  distributions  of 
the  constituent  parts  and  strength  factors  be  known.  It  may,  for  example, 
be  necessary  that  the  distribution  of  the  dimensions  of  depth,  beam,  and 
the  area  of  flanges  be  known.  Such  quantities  are  much  more  amenable 
to  scrutiny  in  ships  than  the  overall  testing  of  the  hull  girder.  As 
discussed  in  Section  7*0,  however,  not  much  data  presently  exist  for  many 
of  the  variables,  and  consequently  the  distributions  themselves  cannot 
be  identified.  This  would  seem  to  be  a  promising  area  in  the  future, 
if  an  effort  is  made  to  collect  such  data. 

If  the  distribution  of  the  constituent  parts  and  functions 
are  known,  there  are  various  methods  for  synthesizing  their  distribution 
to  obtain  the  overall  strength  distribution.  Reference  [7]  gives  eight 
methods: 


0  The  algebra  of  normal  function  method 
0  The  change  of  variable  method 
0  The  moment  generating  function  method 

0  The  Fourier  transform,  convolution,  and  inversion  method 
0  The  Mel  1  in  transform,  convolution  and  inversion  method 
0  The  characteristic  function  method 
0  The  cumulative  distribution  function  method 
0  The  Monte  Carlo  method 

The  Monte  Carlo  method  will  always  give  results  even  for  complex  functions 
of  non- i dent ical ly  distributed  random  variables  although  the  length  and 
complexity  of  the  computations  will  reportedly  be  quite  extensive  and 
possibly  unrealistic. 

The  third  approach  requires  that  assumptions  be  made 
concerning  the  distribution  of  the  strength.  Of  course  the  same  could  be 
done  with  the  constituent  parts  and  factors,  and  the  second  method  used 
to  construct  the  strength  distribution.  This  approach  seems  to  be 
universal  in  the  literature  for  civil  engineering  and  naval  architecture. 

It  is  natural  that  these  two  disciplines  would  make  greater  use  of  this 
last  approach  because  of  the  size  and  complexity  of  the  structure  analyzed. 


-10- 


This  approach  requires  the  adoption  of  a  distribution  (such 
as  the  normal,  lognormal,  Weibull,  etc.)  and  the  specification  of  necessary 
parameters  of  the  distribution  to  obtain  numerical  values  from  tabulated 
density  and  distribution  functions.  The  necessary  parameters  are  at  least  the 
first  and  second  moments  of  the  distribution,  the  mean  and  variance. 

Most  of  the  assumed  distributions  in  the  literature  on 
structural  analysis  are  the  normal  and  the  lognormal  distributions.  It 
would  seem  natural  for  investigators  to  make  such  assumptions  since  expe¬ 
rimental  measurements  in  science  and  engineering  seem  to  approximate, 
rather  well,  the  norma)  law.  However,  the  integrations  of  Equations  (6) 
and  (7)  for  the  probability  of  failure  involve  important  constituent 
parts  at  the  tail  endof  the  distributions  which  can  vary  greatly  depending 
on  the  assumed  distributions.  In  reference  [4],  it  is  stated  that  for 
the  probability  of  failure  10  ,  the  calculated  probability  is  sen¬ 

sitive  to  the  assumed  distribution  and  the  results  can  only  be  used  _ 
relatively.  On  the  other  hand  for  probabilities  of  failure  >  10  , 

such  problems  would  not  be  too  serious. 

As  reported  in  [ 1 5l »  the  record  of  world  ship  catastrophes 
indicate  a  current  probability  of  failure  for  ships  in  the  order  of 
10  ^  so  that  these  approximations  may  not  be  a  problem  in  the  case  of 
ships  if  the  historical  safety  levels  are  considered  adequate. 

If  the  strength  is  assumed  to  be  normally  distributed, 
the  probability  density  and  distribution  functions  are: 


4(s) 


exp  -l/2(^-^) 


*5  ts)  -ffs(sUs  -fs 


where:  m,.  =  mean  of  strength  S 

=  standard  deviation  of  strength  S 

»  standard  tabulated  normal  function. 

Consequently,  under  such  an  assumpt ion, the  only  quantities  that  need  to  be 
determined  are  the  means  and  variances  of  the  strength.  Then,  the  pro¬ 
bability  of  failure  given  by  Equations  (6)  and  (7)  can  be  evaluated 
(provided  the  load  distribution  is  known).  The  latter  statement  is  not 
trivial  since,  in  fact,  the  means  and  variances  of  ships  strength  are  not 
easily  determinable. 

The  approach, in  general,  has  been  to  expand  the  strength 
function  in  terms  of  its  constituents  in  a  Taylor  Series  about  the  means 
of  the  constituents: 

S-f(«y  e2>— -»En) 


.  ,(r.  (-3^.  (22) 

+  1/2 Z  (c.  “  T  )  (X~)  +  . —  +  (Remainder) 

**/  '  1  d~'2£i 

in  which  tne  derivatives  are  evaluateo  at  tne  constituent  means,  r. ,  r  , -- 

— €  and  the  remainder  consists  of  the  higher  derivatives.  4 

n 

If  it  is  assumed  that  the  higher  derivatives  are  small 
or  zero  and  that  the  coefficients  of  variation  of  the  constituents  are 
small,  in  the  order  of  15  per  cent  or  less  [15],  then  Equation  (22)  can 
be  linearized  and  the  following  obtained: 

»s  2  ’  c2  ’  ^  1(23) 


at  ^1 


2  f3f  » 2  _  r 
1  ^  E.  {+j. 
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Where  is  the  correlation  coefficient  between  e.  and  £j. 

These  assumptions  may  Tiot  turn  out  to  be  correct  for  all  ships  for  all 
modes  of  failure.  It  is  indicated  in  [Hi]  that  the  inclusion  of  non- 
linearities  in  the  strength  distribution  causes  various  changes  only  in 
the  predictions  of  long-term  probability  of  failure. 


Further,  making  the  assumption  that  the  constituent  parts 
are  statistically  independent,  the  correlation  becomes  zero  and  Equation  (2*0 
reduces  to: 


77- 
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(25) 


Equations  (23)  and  (25)  have  been  used  in  ship  studies  to 
date.  The  assumption  of  zero  ca-relation  inherent  in  Equation  (25)  may 
be  reasonable  for  many  of  the  constituent  parts.  For  example,  in  the 
case  of  the  strength  defined  by  Equation  (27), the  beam  (B  )  should  have 
no  effect  on  the  depth  (D)  and  similarly  both  D  and  B  should  have  no 
effect  on  plate  thicknesses  tf  and  t  On  the  other  hand,  as  an  example, 
the  strengths  in  different  failure  modes  of  the  same  panel  may  be  highly 
correlated  [56], 

If  Equation  (25)  is  written  in  terms  of  a  coefficient  of 

variation  (COV) : 


where: 


gs  9  strength  COV 


COV's  of  constituent  parts 
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Equations  (23),  (25),  and  (26)  then  give  the  strength  parameter's  mean, 
var  ianca, and  COV  respectively  in  terms  of  the  means  and  variances  of 
the  constituent  parts,  (  e  ).  These  must  be  determined  from  data  or  by 
estimation  as  discussed  in  detail  in  Section  7-0-  The  definition  of 
the  strength  is  then  complete  and  the  probability  of  failure  can  then 
be  evaluated.  The  greatest  amount  of  effort  is  needed  in  determining 
the  strength  COV ,  and  is  only  a  fraction  of  that  required  by  the  first 
two  approaches.  One  would,  of  course, have  a  lesser  degree  of  confidence 
in  the  results. 

The  fourth  approach  requires  only  that  the  COV  or  the 
mean  and  variance  of  the  strength  be  known.  The  procedure  to  obtain 
these  was  just  given  above.  These  data  can  only  be  used  in  the  semi- 
probabilistic  methods  outlined  in  Sections  3.2.3  and  3.2.4.  This  approach 
requires  the  least  computational  effort  to  obtain  its  results. 

3.3.4  Time  Dependent  Strengths 

In  general,  whenever  a  critical  failure  mode  involves 
a  time  variant  strength  such  as  it  does  in  the  cases  of  fracture,  fatigue, 
thermal  effect,  and  corrosion,  the  variations  with  time  must  be  accounted 
for.  If  the  strength  can  be  treated  as  a  function  of  time,  the  general 
probabilistic  procedures  presented  previously  can  be  utilized. 

Mechanical  reliability  for  components  exposed  to  fatigue  is 
discussed  in  [8]  and  [9].  The  approach  therein  is  to  use  the  form  of 
strength  given  by  Equation  (19)  which  would  take  care  of  some  time-dependent 
effects  through  the  K  coefficients;  this  is  implied  but  not  stated. 

From  the  standpoint  of  fatigue,  the  following  problems  ore 
directly  addressed  in  these  references: 

6  Fatigue  under  a  fixed  alternating  load  level,  given 
the  "cycles  to  failure’’  distribution  of  the  component. 

0  Fatigue  for  a  specified  life  given  the  broad  band 
strength  and  load  distributions  for  that  life. 

0  Cumulative  fatigue  under  sequential  groups  of  stresses, 
each  group  having  a  specific  number  of  cycles  and  the 
same  maximum  and  mean  alternating  stress  levels. 

The  approaches  to  solving  these  problems  are  identical  to  those  previously 
discussed  herein  in  that  all  analyses  are  performed  at  a  given  tine  in  the 
life  of  the  component  and  at  a  constant  load  level. 

Reference  [17]  reports  on  studies  conducted  to 
investigate  time-varying  structural  probabilistic  strengths  in  the  jet 
engine  field.  The  basis  of  the  general  procedure  proposed  is  a  compu¬ 
tational  sequence  to  determine  probability  of  failure  vs  time  consisting 
of  two  phases:  the  first  is  a  failure  probability  phase  and  the  second 
a  deqradation  of  strength  phase.  Thus,  a  probability  of  failure  calcu¬ 
lation  is  made,  followed  by  a  strength  degradation  calculation  reflecting 
some  operation  time.  The  sequence  can  be  repeated  indefinitely.  The 
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crux  of  the  procedure  revolves  around  identifying  a  time-varying  strength 
degradation  scenario-  Several  types  are  proposed  but  the  analyses  reported 
in  that  paper  were  of  a  "preliminary"  nature.  It  is  noted  that  additional 
work  was  in  progress  at  that  time. 


During  the  course  of  the  study  presented  herein,  a  po¬ 
tential  scenario  for  corrosion  of  ship  hulls  was  envisioned.  If  the 
mode  of  failure  under  consideration  is  that  of  yielding  during  bending 
of  the  hull  girder  as  a  "free-free"  beam,  it  can  easily  be  shown  that 
the  strength  equation  is: 
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A^/2B  =  equivalent  thickness  of  one  flange 

A  /2D=  equivalent  thickness  of  one  web 
w 


if  corrosion  is  introduced,  then  A^.,  A^,  t^.  and  t 
become  functions  of  time  as  the  plating  corrodes.  "v  w 


The  plate  thicknesses  may  then  be  considered  a  function  of  time  as  follows: 

t(t)  -  tQ  -rMT>  (28) 

where:  t(x)  =,  Plate  thickness  in  time  ~ 


tD  =*  Original  thickness  at  r=>o 

Rc  **  Rate  of  corrosion,  also  a  random  variable 

The  strength  would  then  become  a  function  of  time  and  the  probability  of 
failure  could  be  estimated  at  various  times  during  the  ship's  life  using 
the  probabilistic  theory  previously  presented.  A1 ternat i vely, the  original 
strength  at  time  t=o  could  be  multiplied  by  a  factor  kc, ref iect ing  equation 
(?3) ,  also  a  random  variable,  to  account  for  a  specific  reduction  in 
strength  at  a  certain  time  in  the  vessel  life. 

Another  approach  to  consider  the  effect  of  corrosion 
which  does  not  result  in  a  time  dependent  strength  is  to  take  the  total 
plate  thickness  as  the  sum  of  the  thickness  required  for  limiting  stresses. 
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t  ,  plus  a  thickness  for  corrosion  allowance,  t  [62]: 
n  c 


Which  by  Equation  (26)  yields: 

62  =  (l£_)2  •  6?  +  (I^-)2  62  (30) 

t  t  n  1  c 

v.’here : 

&  =  COV  of  the  plate  thickness  due  to  production  tolerances 

*n 

6  =  COV  due  to  corrosion 

c 

As  pointed  out  in  [62],  the  corrosion  rate  will  vary  from  one  group  of 
strength  members  to  another  and  this  has  been  addressed  by  others  using 
a  Monte  Carlo  simulation  technique  [65]. 

In  [19]»a  method  is  presented  for  probabilistic  analysis 
of  fatigue-crack  initiation  at  a  butt-welded  joint.  The  procedure  is  used 
for  analyzing  both  the  longitudinal  and  transverse  structural  members  of 
a  tanker  subjected  to  random  still  water  and  wave  loads.  This  reference 
represents  the  only  source  found  during  the  course  of  this  study  which 
gives  a  probabilistic  evaluation  of  ship  structure  fatigue.  The  strength 
function  given  therein  is  based  on  Miners'  law  and  on  the  coefficients 
of  a  logarithmic  linear  approximation  of  the  S-N  curve,  which  are  regarded 
as  random  variables.  A  sensitivity  analysis  on  these  random  variables 
is  also  presented.  The  degradation  of  strength  in  time  by  factors  other 
then  fatigue  is  not  considered  and  it  is  noted  that: 

"because  of  lack  of  sufficient  amount  of  statistic  data 
or  quantitative  information  on  unexpected  defects  in  hull  structure,  this 
study  is  limited  to  within  a  range  of  treating  only  a  standard  ship  which 
is  built  through  sound  workmanship  of  well  quality-controlled  fabrication 
and  is  put  into  service  with  satisfactory  maintenance  under  normal  ope¬ 
rating  conditions.  It  should, therefore, be  clearly  born  in  mind  that 
the  results  obtained  by  this  analysis  will  provide  information  on 
the  reliability  of  ship  structures  merely  on  the  basis  of  design-oriented 
point  of  view."  [19]. 

The  approach  used  in  [5],  as  previously  discussed  in 
Section  3-2.4.,  has  been  extended  therein  to  constant  stress  range  and 
random  fatigue. 

3-4  Load  Statistics 
3.4.1  Genera  1 


As  discussed  in  Section  2.0,  the  objectives  of  this  study 
do  not  include  details  concerning  the  load  distribution.  However,  since 
the  load  is  one  of  the  two  major  considerations  of  any  probabilistic 
structural  design,  it  will  be  discussed  here  from  the  standpoint  of 
characteristics  and  mechanics  that  must  be  considered  for  application  in 
probabilistic  structural  design.  The  literature  on  loads  does  not 
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address  this  point  extensively.  A  qualitative  appraisal  of  the  situation 
with  respect  to  loads  applied  to  ships  is  included  in  Section  5-0  of 
this  report. 


The  types  of  loads  applied  to  the  hull  girder  consist 
or  the  following  (59l  ; 

°  Calm  water  due  to  weight  and  buoyancy. 

°  Ship's  own  wave  train. 

°  Thermal  effects. 

°  Quasi-static  wave  induced  (low  frequency). 

°  Dynamic  (high  frequency):  including  slamming, 

whipping,  springing,  and  propeller  induced  vibration. 

3.4.2  Equations  and  Distributions 

Equations  (1)  through  (7)  deal  with  expressions  for 
the  probability  of  failure,  reliability,  and  margin  of  safety.  In  these 
expressions  strength  and  load  carry  the  same  weight  and  require  the  same 
type  of  expressions  for  their  mathematical  description.  Hence,  all  that 
has  been  stated  for  the  strength  equations  and  distributions  would  apply 
in  most  cases  to  the  load  distributions  as  well. 

With  respect  to  ships,  the  procedures  of  synthesizing 
distributions  of  the  constituent  parts  into  that  of  the  whole  should  be 
emphasized.  The  procedure  for  combining  still  water  and  wave  bending 
moments,  springing,  slamming,  and  thermal  effects  should  be  similar  to 
that  presented  in  3-3-3  for  strength  distributions. 

The  analyses  to  be  found  in  the  literature  on  probabilistic 
structural  design  of  ships  have  only  considered  still  water  and  wave 
bending  moments  directly.  This  is  primarily  due  to  lack  of  information 
applicable  to  other  types  of  loads,  as  discussed  further  in  Section  5-0. 

It  should  be  pointed  out  here,  however,  that  in  any  complete  probabilistic 
analysis,  the  total  load  must  be  considered. 

In  the  case  of  longitudinal  strength,  this  total  load  will 
include  the  effects  of  local  loadings,  such  as  that  due  to  water  head, 
since  this  will  add  a  random  load  toward  increasing  the.  overall  load  and 
hence,  the  stress. 


With  respect  to  specific  distributions  proposed  in  the 
literature,  those  found  in  [10]  have  been  used  in  probabilistic  structural 
analyses  of  ships  presented  therein;  the  wave  bending  moments  and  still- 
water  bending  moments  have  been  considered.  The  amplitudes  of  the  wave 
bending  moments  are  assumed  to  follow  a  Rayleigh  distribution  in  the  short 
term,  and  an  exponential  probability  law  in  the  long  term.  Using  the 
Weibull  distribution,  both  the  short-term  and  the  long-term  wave  distri¬ 
bution  and  density  functions,  respectively,  are  given  as  follows: 


f|>> 

Fl  (x) 
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(30 

(32) 


where:  l  m  2  for  short  term 

“  1  for  long  term 
k  ■  /E~ for  short  term 

k  **  A  for  long  term 

E  ■  mean  square  value  of  L  taken  over  a  short  period  of  time 
A  =  expected  value  of  L  taken  over  a  long  period  of  time 

It  should  be  pointed  out  that  in  [60]  it  is  shown  that  the 
exponential  law  underestimates  the  data  measured  onboard  an 
Ore/Bulk/Oil  carrier.  Therein, it  is  concluded  that  mathematical  models 
based  on  the  normal  or  general  Weibull  distributions  give  excellent  agree¬ 
ment  with  statistical  data  for  the  ship  analyzed.  Reference  [6l]  shows 
that  for  two  other  ships*  the  Weibull  distribution  does  not  exactly  fit 
the  data. 


a 


In  Reference  [12],  "order  statistics"  are  used  to  obtain 
the  extreme  wave  bending  moment  density  and  distribution  functions  using 
equations  (31)  and  (32).  These  extreme  functions  become: 

*(y)  „  Hi  (y/k) 1  - e-  <y/k> £ 1 1  -e" (y/k) n'1  y>°  (33) 

Yn  k 

*(y)  =-  P[Y«y]  =  Il-e“(y/k)t]n  y*0  (34) 

Yn  n 

where  n  is  the  number  of  wave  records  considered. 


The  still-water  bending  moment  is  incorporated  first  as 
deterministic  and  then  as  a  normally  distributed  random  variable.  The 
combined  still-water  and  wave  bending  moment  probability  density  and 
distribution  functions  in  the  deterministic  case  are: 
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-  0,  otherwise, 
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0- otherwise,  « 

The  probsbil ity  density  and  distribution  functions  in  the 


where  m  is  the  deterministic  bending  moment, 
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and 
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where  m  and  o  are  the  mean  and  standard  deviation  of  the  still  water  bending 
moment  respectively. 
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SECTION  4.0 


MOOES  OF  HULL  FAILURE 


4. 1  General 

It  is  well  known  that  the  design  of  a  ship's  hull  girder 
from  the  standpoint  of  longitudinal  strength  is  usually  performed  by 
considering  yield  failure  of  the  hull  girder  as  a  free-free  beam  in  bending. 

The  load  is  normally  determined  by  balancing  the  ship  on  an  "extreme 
wave"  for  both  hogging  and  sagging  conditions  and  the  resulting  stress 
must  remain  below  an  allowable  level.  Factors  of  safety  based  on 
experience  are  contained  in  the  loads  and  the  allowable  stresses.  Ex¬ 
perience  has  shown  that  such  an  approach  leads  to  probabilities  of 
commercial  ship  failures  in  the  order  of  10  3,  [15,  5 9 1  »  although  the  modes  of 
the  failures  are  not  all  known. 

In  turning  to  probabilistic  structural  design,  as  pointed 
out  in  Section  3-0»  all  conventional  factors  of  safety  must  be  stripped  away 
and  accurate  distributions  of  load  and  strength  must  be  determined. 

Further,  all  potential  modes  of  failure  must  be  analyzed  in  separate 
calculations. 


This  last  aspect  may  appear  subtle  to  some;  but  one  must 
remember  that  the  yield  failure  of  the  hull  girder  as  a  beam  is  not 
the  only  potential  mode  of  failure  of  a  ship  hull  girder.  V/ith  the 
historical  conventional  factor  of  safety  approach  on  this  yield  failure 
mode,  other  modes  of  failure  may  also  be  automatically  taken  care  of 
but  with  smaller  margin  and, therefore, wi th  less  of  an  effective  factor 
of  safety.  This,of  course, is  the  major  shortcoming  of  the  conventional 
factor  of  safety  approach  and  is  rooted  in  its  empiricism. 

Consequently,  In  the  probabilistic  structural  analysis,  all 
the  potential  modes  of  failure  of  the  hull  girder  must  be  identified  and 
analyzed.  The  output  may  again  be  a  factor  of  safety,  but  its  determination 
would  be  on  a  more  rational  basis. 

4 . 2  Modes  of  Fa i lure  of  the  Hul 1  G i r de r 


Modes  of  failure  of  the  huli  girder  from  a  longitudinal  strength 
standpoint  can  be  grouped  into  the  following: 

*  Yield  failure  due  to  bending  of  the  ship  considered 
as  a  free-free  beam 

°  Compression  instability  buckling 
°  Brittle  fracture 

0  Fatigue  fracture 

#  Ultimate  plastic  collapse 
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As  previously  stated,  longitudinal  strength  in  hull  girder  design 
is  usually  based  on  the  deterministic  evaluation  of  beam  bending  with  factors  of 
safety  to  prevent  a  yield  failure.  However,  it  is  interesting  to  note  that 
various  investigators  have  indicated  this  not  to  be  the  most  significant 
mode  of  failure,  [21]  and  [23]. 

In  [21],  it  is  shown  that  compressive  and  tensile  strengths  of 
even  poorly  built  ships  are  adequate  to  withstand  the  most  severe  wave 
bending  moments.  With  respect  to  brittle  fracture,  it  is  noted  that 
fractures  cannot  initiate  because  the  quality  of  workmanship  today  is 
high  and  the  nominal  stresses  are  usually  low.  However,  if  higher  allowable 
stresses  in  hull  materials  are  used,  then  means  of  arresting  cracks  will 
have  to  be  considered.  The  feeling  is  that  the  brittle-fracture  problem 
can  be  eliminated  by  proper  use  of  crack-arresting  steels  under  any 
circumstances.  In  the  future,  the  problem  may  be  restricted  to  fatigue 
cracks  and  how  large  they  may  be  allowed  to  get  without  leading  to 
unstable  fracture.  Fracture-mechanics  investigations  are  proposed 
for  this  analysis.  A  statement  made  in  (21]  is  of  interest: 

"So  much  for  the  brittle-fracture  problem.  It  is  quite 
possible  that  within  10  or  20  years  it  has  disappeared  from  shipbuilding. 

Then  the  level  of  permissible  stresses  will  be  to  a  large  extent 
determined  by  fatigue  considerations.  In  fact  it  docs  so  already  now¬ 
adays  together  with  brittle  fracture,  buckling  of  bulkheads  and  webs 
of  deep  frames  and  bottom  damage  due  to  slamming.  It  seems  that 
not  everyone  is  aware  of  this  fact.  There  are  even  investigators, 
dedicating  their  time  to  wave  bending  moments,  who  are  not  much 
interested  in  fatigue." 

In  [22],  a  method  is  presented  for  the  determination  of  the 
ultimate  plastic  moment  of  the  hull  girder.  It  is  stated  that  elastic 
stresses  from  the  conventional  approach: 

"may  be  influenced  by  residual  reaction  or  thermal  effects 
to  such  an  uncertain  extent  that  the  stresses  thus  calculated  are  some¬ 
times  regarded  as  having  only  comparative  rather  than  absolute  value. 

The  ultimate  strength  of  a  ship  is  likely  to  be  influenced  by  these 
uncertainties  to  a  much  smaller  extent,  so  that  the  calculated  hull 
bending  moment  should  give  a  reliable  indication  of  thq  true  bending 
strength  of  the  hull.  It  must  be  emphasized  again,  however,  that  the 
possibility  of  premature  failure  by  major  hull  fracture  must  be  guarded 
against  by  proper  design  and  construction  details  and  control  of  material 
quality.  If  this  is  true,  then  overall  hull  girder  failure  can  only 
occur  through  yielding  and  buckling,  in  the  way  assumed  in  this  analysis." 

However,  discussions  of  the  cited  reference  indicate  that 
buckling  has  been  eliminated  to  a  very  great  extent  and  brittle  fracture 
is  the  principal  hazard,  [23],  and  that  low-cycle  fatigue  leading  to 
local  failure  and  hastening  the  complete  "breaking  its  back"  before 
the  ideal  ultimate  failure  load  is  the  primary  problem,  [2*0  . 

It  is  proposed  in  [11]  that  the  fracture  modes  of  failure 
can  be  avoided  providing  care  is  taken  in  material  selection  and  inspec¬ 
tions  are  made  periodically.  In  conclusion,  it  is  in  effect  stated  that 
only  adequate  safeguard  against  the  occurence  of  plastic  collapse  need 
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be  provided.  This  is  tantamount  to  the  considerations  of  compression 
instability  which  are  summarized  as  strut-panel  and  tripping  of 
stiffeners  locally  as  well  as  overall  grillage  buckling. 

In  [18],  the  importance  of  analyzing  various  modes  of 
failures  and  damage  to  ship  structures  is  pointed  out.  Results  of  a 
20-year-lifetime  probability  analysis  are  given  with  respect  to  yielding, 
local  buckling,  total  collapse,  and  fatigue  type  failures.  Effects  of 
local  water  pressure  are  also  included.  The  results,  quoting  from 
the  afore-mentioned  study,  were  that  uThe  probability  of  fatigue  crack 
initiation  is  comparatively  high,  whereas  for  ductile  failures, 
probability  of  local  collapse  of  bottom  longitudinals  is  fairly  significant, 
followed  by  the  yield  failure  of  deck  or  bottom  plating,  and  very  low 
probability  of  total  plastic  collapse  of  the  hull  girder." 

As  a  further  complication  to  the  problem,  one  must  also 
remember  that  many  of  the  proposed  modes  of  failure  have  been  inves¬ 
tigated  from  a  "stress  at  a  point"  view  and  due  to  primary  hull  stresses 
only.  However,  the  hull  girder  has  the  capability  of  redistributing 
stresses  once  it  yields  at  a  point.  The  total  principal  stress  must, 
therefore, be  determined  by  the  superposition  of  primary,  secondary, 
and  tertiary  stresses.  Again  these  considerations  are  not  important 
in  the  usual  empirical  approach  to  longitudinal  strength  but  are  of 
great  concern  in  any  precise  structural  analysis. 

i» .  3  Conclusion 


It  is  obvious  from  the  foregoing  that  the  mode  of  failure 
for  a  ship  hull  girder  is  not  specifically  known.  In  fact,  it  seems 
perfectly  plausible  that  the  mode  of  failure  may  vary  depending  on 
the  design  as  is  generally  experienced  in  structural  design.  Furthermore, 
the  overall  probability  of  failure  requires  that  all  probabilities 
of  failure  of  individual  modes  of  failure  be  known  and  combined,  and 
that  the  total  stresses  including  any  local  stresses  must  be  considered. 

As  opposed  to  this  situation,  in  the  examples  of 
probabilistic  structural  design  for  machine  parts,  such  as  the  one 
In  [8],  the  mode  of  failure  and  various  stress  components  acting  on 
the  parts  are  exactly  known;  and  it  is  emphasized  that  fhis  must  be 
che  case.  Section  3-0  discusses  this  subject  in  more  detail. 
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SECTION  5-0 


LOADINGS 


The  present  study  is  not  concerned  with  any  investigation 
of  loadings  on  the  hull  girder  other  than  to  obtain,  from  a  review  of 
the  literature,  data  on  loadings  needed  to  perform  an  example  calculation. 
Yet,  this  point  must  be  addressed  in  principle,  since  it  shares  an 
equal  portion  with  the  strength  of  structures  in  the  orobabi 1 i st i c 
structural  design  theory.  In  other  words,  in  order  to  perform  pro¬ 
babilistic  structural  analyses,  all  must  be  known  about  both  the  load 
and  the  strength. 

It  was  stated  in  Section  3>0  that  the  load  considered  must 
be  the  total  load  acting  on  the  structure  to  cause  the  particular 
failure  in  question.  In  relation  to  the  longitudinal  strength  of 
a  ship  hull  girder,  such  loads  would  include  still-water  bending  moments, 
wave-induced  bending  moments,  springing  induced  bending  moments, 
slamming  induced  bending  moments  of  all  types,  transient  deck  loads 
due  to  weather,  thermal  effects,  and  bending  moments  due  to  the  ships 
own  wave  train.  Except  for  the  wave  loads,  there  is  very  little 
In  the  literature  concerning  the  statistical  data  for  these  various 
loadings. 


From  the  standpoint  of  analytically  determining  lifetime 
wave  loads.  Reference  [12]  presents  a  procedure  for  determining  the 
extreme  values  of  the  wave  bending  moment  using  "order  statistics" 
and  assuming  that  the  distribution  of  the  maximums  is  of  the  Weibull  type. 

Several  investigators  have  presented  statistical  full-scale 
data  measurements  of  wave  bending  moments  for  actual  vessels,  [25],  [26], 
[27].  In  the  measurements  presented,  the  effects  of  springing  and  whipping 
were  filtered  out.  The  results  are  curves  of  cumulative  long-term 
distribution  of  the  average  bending  moments  whicn  show  the  probabilities, 
per  cycle  of  load  [27],  of  exceeding  different  levels  of  these  bending 
moments  during  a  ship's  lifetime.  Figure  1  is  an  example  reproduced 
from  [26].  A  method  for  converting  these  loads  per  cycle  to  a  cumulative 
probability  curve  for  the  ship's  lifetime  is  indicated  in  [28].  Following 
this  procedure,  a  form  of  long  term  distribution  must  be  assumed. 

As  discussed  in  Section  3-**,  the  different  assumed  long-term-distribution 
shapes  fit  the  measured  data  differently  [10,  60,  61]. 

There  does  not  appear  to  be  enough  data  nor  any  analytical 
methods  in  the  literature  for  determining  the  statistical  distribution 
of  the  other  loads  mentioned  above.  The  ship  structural  reliability 
studies  presented  in  [10]  assumed  both  deterministic  and  normally 
distributed  still-water  bending  moments.  Reference  [29]  discusses  the 
computations  of  wave  slamming  and  springing  bending  moments  in  the 
context  of  a  probabilistic  structural  analysis,  but  it  is  pointed  out  that 
much  verification  must  be  made  with  respect  to  slamming  and  springing 
before  the  procedures  can  be  used.  It  is  also  noted  that  with  regard 
to  the  structural  probabilistic  analysis,  springing  and  slamming  were 
not  incorporated  although  they  might  easily  be. 
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SECTION  6.0 


PROBABILISTIC  STRUCTURAL  ANALYSIS  OF 


6.1  General 


SHIP  HULL  LONGITUDINAL  STRENGTH 


All  of  the  most  essential  considerations  for  probabilistic  structural 
design,  discussed  in  preceding  sections,  would  apply  to  transverse  and 
torsional  hull  strength  as  well  as  longitudinal.  However,  transverse  and 
torsional  hull  strength  are  beyond  the  scope  of  this  study  (Section  2.0). 

It  is  clear  that  there  are  significant  problems  concerning  the 
data,  theory,  and  techniques  that  stand  in  the  way  of  a  completely  rational 
probabilistic  hull  girder  longitudinal  strength  analysis.  This  is  to  say 
that  the  probabilities  of  failure  from  such  an  analysis  could  only  be  used 
In  a  relative  sense;  and  even  then  the  comparison  of  modes  of  failure 
might  be  questionable  due  to  possible  better  input  to  one  mode  of  failure 
analysis  than  the  other. 

Other  investigators  have  discussed  this  point.  It  is  stated  in 
[18]  that  the  relative  assessment  of  probabilities  of  failure  may  be  one  of 
the  useful  methods  of  evaluation  of  ship  structures.  In  Reference  [13],  one 
of  the  motivations  behind  the  approximate  approach  presented  therein  was 
that  "probabilistic  analysis  of  structural  safety  for  ships  is  difficult  at 
the  present  time  because  the  available  data  are  too  limited  to  provide  the 
exact  forms  of  the  probability  distributions  of  the  bending  moment  and 
the  ship  strength."  Reportedly,  the  sample  size  required  is  of  the  order 
of  multimillion  pieces  of  records  or  data  [30].  Two  more  recent  papers, 

162]  and  [63],  also  discuss  this  point. 

One  other  aspect  of  probabilistic  design  which  has  received  mention 
but  net  much  analysis  is  the  problem  of  determining  the  acceptable  limiting 
value  of  the  probability  of  failure.  It  was  mentioned  previously  that  the 
current  level,  based  on  actual  occurrences,  was  determined  in  a  study  [15,  59]- 

The  two  emerging  problems,  l.e.  the  lack  of  available  data  and 
techniques  to  perform  an  accurate  probability  of  failure  analysis  and  the 
absence  of  an  acceptable  limit  to  the  probability  of  failure,  point  to  a 
need  for  the  following  three  overall  efforts: 

°  Continue  to  develop  techniques  and  obtain  data  for  both 
load  and  strength  for  probabilistic  analysis. 

°  Perform  absolute  probability  of  failure  analyses  for 
different  ships,  compare  and  update  the  results  as 
better  data  and  techniques  are  developed. 

°  From  the  data  presently  available  on  ship  failures  of  all 
types  for  all  types  of  ships,  perform  semi -probabi 1 i st i c 
analyses  to  identify  safety  factors  of  current  and 
past  ships. 

The  first  of  the  above  is  needed  for  the  advancement  of  probabilistic 
structural  analysis  methodology. 
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The  second  is  necessary  to  obtain  results  using  the  latest  collected 
data  and  the  theoretical  and  testing  methods  in  a  developing  analysis 
procedure,  and  to  compare  these  results  with  what  is  known  and  with  each 
other.  Ultimately,  the  results  of  such  a  procedure  could  be  used,  in 
conjunction  with  nun-structural  aspects,  to  obtain  an  acceptable  limit  for 
the  probability  of  failure.  This  would  lead  to  the  determination  of  the 
resulting  factors  of  safety  to  be  used  by  designers. 

The  third  arid  last  effort  is  needed  to  reap  immediate  fruits  from 
the  probabilistic  design  approach.  As  mentioned  earlier,  current  longitudinal 
strength  procedures  only  consider  ductile  yielding  of  the  hull  girder  due 
to  an  equivalent  wave  imposing  a  vertical  bending  moment.  An  extensive 
analysis  of  ships,  particularly  those  that  have  failed  in  longitudinal 
strength,  considering  all  the  modes  of  failure  and  using  known  loads,  can 
produce  a  better  understanding  of  which  modes  of  failure  are  most  significant, 
what  the  factors  of  safety  are  for  these  modes,  and  possibly  indicate  trends 
with  respect  to  ship  type,  size,  area  of  operation,  etc. 

It  should  be  noted  that  the  factor  of  safety  discussed  in  conjunction 
with  the  second  effort  above  is  different  from  that  discussed  in  the  third. 
Factors  of  safety  that  come  from  an  exact  probabilistic  analysis  are  to  be 
based  on  the  exact  knowledge  of  the  load  and  strength.  This  would  be  the 
case  with  the  second  item  but  in  the  case  of  the  third,  uncertainty  concerning 
the  load  and  strength  would  be  tied  into  that  safety  factor  but  it  should  involve 
less  uncertainty  than  in  current  procedures.  As  more  becomes  known  about 
the  strength  and  load,  the  results  of  the  third  item  can  of  course  be  updated. 

The  proposed  approaches  for  the  above  three  areas  of  effort  are 
presented  below  In  greater  detail. 

6.2  Development  of  a  Probabilistic  Structural  Analysis  Methodology 

In  Sections  3.  **  and  5,  the  areas  were  identified  where  more 
theoretical  studies  need  to  be  performed  and  data  collected  from  full-scale 
experimental  measurements.  These  areas  can  further  be  divided  into  strength 
and  load  distributions,  strength  equations,  and  time-dependent  analyses. 

6.2.1  Strength  and  Load  Distributions 

As  discussed  In  Section  3-0,  if  the  exact  form,  and  the 
magnitudes  of  the  strength  and  load  distributions  are  known,  it  is  simple 
to  determine  the  probability  of  failure.  The  problem  is  that  for  both 
strength  and  load,  these  distributions  do  not  exist.  The  problems  for  load 
•re  discussed  in  Section  5-0,  and  will  not  be  elaborated  upon  here  since  this 
Is  not  the  specific  area  addressed  by  this  project  and  is  being  addressed 
elsewhere  [32], 


In  the  case  of  the  strength,  the  procedure  has  been  to 
synthesize  the  distribution  either  by  estimating  only  the  coefficients  of 
variation  of  strength  variables,  l.e.  constituent  parts  whose  distributions 
are  not  known,  or  by  making  an  assumption  as  to  what  type  of  distribution  the 
strength  follows.  Higher  level  syntheses  discussed  in  Section  333  consist 
of  the  determination  of  strength  statistics  from  the  assumed  distribution  of  the 
constituent  parts,  from  the  known  distribution  of  constituent  pans,  or  the 
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determination  of  strength  statistics  from  actual  testing.  The  latter  would 
give  the  greatest  accuracy  but,  as  discussed  previously,  is  very  nearly 
impossible  due  to  its  extent. 


6.2.2  Strength  Equations 

A  simple  strength  equation  for  ductile  yielding  due  to  bending 
of  the  hull  girder  is  given  in  Section  3-0,  Equation  (27)-  The  whole  structural 
description  rests  with  equations  like  these  and  the  variables  they  include, 
since  it  is  the  statistics  of  these  variables  which  are  used  tc  synthesize  the 
overall  strength  distribution.  Such  an  equation  can  be  lacking  in  the  simplifying 
assumptions  associated  with  its  derivation  or  in  the  number  of  /ariables  it 
contains.  Equations  can  have  the  same  simplifying  assumptions,  which  contribute 
to  subjective  uncertainties  as  discussed  in  Section  7,  but  a  di  ferent  number  of 
variables.  For  instance  consider  the  following: 


S  -  f(M,Sy)  -  Nsy  (39) 

M-  g(0,  tf  ,  B,  tw)  (1»0) 

S  -  f[g(3,  tf  ,  B,  tw),  sy]  *.  D (t pB+1  / 3twD) sy  (<•') 


where:  N  -  deck  or  bottom  section  modulus 

Sy-  tensile  strength 
D  -  section  depth 
8  -  section  beam 
Af-  area  of  flanges 
A*-  area  of  webs 

tf-  Af/2B  -  equivalent  thickness  of  one  flange 
rw-  AW/2D  -  equivalent  thickness  of  one  web 
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Equation  (39)  oniy  considers  the  section  modulus  and  tensile  strength  and  is 
based  on  the  simplifying  assumptions  of  engineering  beam  theory.  To 
synthesize  a  distribution  from  these,  one  would  need  the  statistics  of  the 
section  modulus  and  of  the  tensile  strength.  The  former  would  be  nearly 
impossible  co  obtain  accurately  since  ships  or  structural  models  would 
have  to  be  tested.  By  breaking  the  section  modulus  into  variables  as  in 
Equation  (40) ,  the  strength  also  becomes  composed  of  more  variables 
amenable  to  more  direct  scrutiny  as  far  as  their  statistics  are  concerned, 
Equation  (4 1 ) .  This  can  naturally  be  extended  to  much  more  subtle 
var ! 3bles . 


The  objective  of  such  an  approach  is  to  define  all  the  variables 
which  can  more  easily  be  measured  and  for  which  statistics  can  be  determined. 
The  methods  discussed  in  Section  3-3.3  should  then  enable  one  to  synthesize 
an  accurate  strength  distribution. 

The  example  of  Equation  (41)  is  usually  intended  for  elastic 
bending  of  the  hull  girder  as  a  beam.  In  the  cases  of  fracture,  fatigue, 
and  ultimate  strength,  more  work  will  probably  have  to  be  done  to  obtain  an 
accurate  expression  for  hull  strength  for  these  modes  of  failure. 

6.2-3  Time -Dependent  Strength  Analyses 

It  was  discussed  in  Section  3.3-4  that  there  is  not  much 
available  in  the  literature  regarding  time-dependent  strengths  of  ships. 

It  i s, therefore, suggested  that  research  be  performed  to  develop  an  accurate 
probability  of  failure  procedure  for  the  analysis  of  time-variant  ship 
strengths,  i.e.  in  the  case  of  fatigue,  corrosion,  etc. 

6. 3  Application  of  Probabilistic  Structural  Analysis  Methodology 

6.3.1  General 

In  order  to  obtain  an  immediate  and  clear  idea  of  ship  longitu¬ 
dinal  strength  from  the  practical  standpoints  of  modes  of  failure,  safety 
margins,  and  probability  of  failure,  the  probabilistic  structural  design 
approach  can  be  used  most  fruitfully  in  a  semi -probabi 1 i st ic  type  of  analysis. 

The  semi -probab i 1 i s t i c  analysis  approach  would  be  the  easiest 
to  apply  and  be  consistent  with  the  assumptions  that  would  be  necessary. 
Applying  a  more  rigorous  and  time-consuming  technique  with  an  equal  amount 
of  additional  assumptions  may  not  add  insight  or  accuracy,  and  it  may  even 
detract  from  the  efforts. 

A  potential  basis  for  such  an  approach  has  been  prfesented  in 
[>3).  It  is  implied  there  that  the  procedure  might  actually  be  used  in 
structural  design.  From  the  standpoint  of  a  designer,  this  seems  highly 
unlikely  in  the  near  future  since  most  ship  longitudinal  strength 
determinations  are  presently  based  on  classification  society  rules,  or 
specifications  which  have  been  developed  through  years  of  experience, 
and  in  all  likelihood  they  would  not  be  changed  until  a  diff&rent  approach 
that  offers  advantage  could  be  established  and  well  proven. 
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However,  in  the  interest  of  obtaining  improved  analysis  methods 
and  insight,  it  is  important  for  researchers  and  designers  to  have  a 
better  appreciation  for  the  basic  structural  phenomena  associated  with 
hull  girder  longitudinal  strength.  As  pointed  out,  the  present  procedure 
of  designing  for  ductile  yielding  of  the  hull  girder  in  vertical  bending 
may  be  mythical  in  some  cases;  and  only  coincidentally,  an  adequate 
strength  for  other  modes  of  failure  may  have  been  accounted  for. 

6.3.2  Method  of  Approach 


The  "Approximate  Probabilistic  Method"  of  Reference  [13]  does 
not  require  assumptions  concerning  the  types  of  distribution  for  load 
and  strength.  This  method  is  a  candidate  for  the  semi -probabi 1 i st i c 
probability  studies  discussed  for  the  short  term.  The  details  of  this 
procedure  are  presented  in  Section  3.2.3- 

The  procedure  was  presented  specifically  for  vertical  bending 
of  the  hull  girder  as  a  beam.  But  it  seems  plausible  that  the  approach 
might  be  extended,  with  some  modification,  to  other  modes  of  failure. 

For  the  specific  purposes  of  obtaining  better  knowledge  on  more  rational 
factors  of  safety  and  on  critical  modes  of  failure,  this  procedure  yields 
the  safety  index  as  given  by  Equation  (11)  of  Section  3-0.  If  this 
index  can  be  evaluated  for  existing  ships,  including  those  that  have 
failed,  for  various  modes  of  failure,  then  the  safety  factors  of  ship 
longitudinal  strength  will  be  better  understood.  The  results  of  limited 
analyses  of  this  type  are  presented  in  [13]  and  [31]  for  oil  tankers  and 
[63]  for  naval  designs.  It  is  of  interest  to  note  that  Figure  2, 
reproduced  from  [13],  gives  an  indication  of  the  probability  of  failure 
for  a  given  safety  index  when  different  distributions  are  assumed. 

The  procedure  for  arriving  at  the  safety  index  would  then 
be  as  fol lows : 

0  Determine  the  mean  strength  for  the  mode  of  failure 
In  question 


ms  -  oNS 

where : 

=  average  failure  stress 

S  =  hull  strength  in  question 
mean  of  strength 

Determine  central  safety  factor. 


0  = 


m2 


where: 

m.,  =  mean  of  total  load 


(42) 


(43) 
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SAFETY  INDEX  Y 


FIG.  2  -  PROBABILITY  OF  FAILURE  VERSUS  THE  SAFETY  INDEX 


°  Determine  safety  index. 

m 


where:  Vs  **  COV  of  strength 
“  COV  of  load. 

The  quantities  tfy  and  S  should  be  obtainable  from  the  design 
calculations  of  the  ship.  m^  can  be  determined  by  measured  data,  empirical 
procedures,  or  by  an  analytical  approach  such  as  the  one  in[;10].  A  correction 
for  loads  not  considered  by  these  approaches  and  the  COV's  of  load 
and  strength  can  be  determined  from  existing  data.  An  example  of  this 
procedure  along  with  a  computer  program  are  presented  in  Section  8.0. 

The  safety  index  Y  will  be  directly  indicative  of  the 
safety  factor,  and  since  it  will  be  derived  from  statistics  of  both 
the  load  and  the  strength,  for  various  modes  of  failure,  it  will  give 
much  more  information  than  the  singular  approach  of  current  design 
based  on  vertical  bending  only. 


SECTION  7.0 


UNCERTAINTIES  IN  HULL  STRENGTH 

7 • 1  Genera  i 


As  mentioned  repeatedly  in  the  preceding  sections,  a  suitable 
approach  to  ship  probabilistic  structural  anal/sis  requires  the  determination 
of  strength  variable  means  and  variances.  The  ratio  of  standard  deviation  to 
mean  is  termed  the  coefficient  of  variation  or,  COV. 

The  statistical  nature  of  the  variables  has  led  to  their  being 
termed  "uncertainties." 

Reference  [ 3 ]  classifies  the  uncerta int ies  into  two  types: 

o  Objective  -  "Measurable  or  quantifiable,  such  as  observed 
statistical  variabilities  and  deductive  probabilistic 
information." 

o  Subjective  -  "No  factual  information  is  available  or  the 
uncertainty  is  not  amenable  to  quantitative  description 
and  must  be  described  and  handled  subjectively  on  the 
basis  of  judgement  and  intuition." 

Once  the  means  and  variances  of  the  objective  and  subjective 
uncerta in,t  ies  have  been  established,  they  can  be  combined  as  follows: 


where:  Vx  =  COV  of  x 

Vxo~  COV  of  objective  uncertainties  of  x 
Vxs*  COV  of  subjective  uncertainties  of  x 

7.2  Objective  Uncertainties 

7.2.1  Genera  1 


The  objective  uncertainties  of  the  hull  longitudinal 
strength  a-e  divided  into  three  grc'ps:  mill  practice,  shipyard  practice, 
and  operat.onal  occurrences. 

o  Mill  Practice 

-  Variation  in  physical  properties  of  materials  including 
ductile,  fatigue  and  fracture  character i s t i cs . 

-  Variation  in  material  thickness  and  shape  dimensions. 

o  Shipyard  Practice 

-  Variations  in  material  scantlings. 

-  Variations  in  fabrication  tolerances. 

-  Variation  in  weldments. 

-  Residual  stresses. 

-Initial  Deflections 
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o  Operational  Occurrences 


-  Corrosion 

-  Wear  and  Tear 

With  the  goal  of  collecting  data  for  these  various 
uncer ta in t ies ,  requests  Were  made  to  Steel  mills  and  the  literature  was 
surveyed.  From  the  steel  mills,  statistical  data  in  the  form  of  either 
distributions,  means  and  variances,  raw  data  points,  or  tolerances  for  at 
least  the  following  quantities  were  requested: 

o  Tensile,  Yield,  and  Ultimate  Strength 
o  Young's  Modulus,  Tangent  Modulus,  and  Shear  Modulus 
o  Duct ili ty 
o  Corrosion  Resistance 

o  Dimensions  of  Manufactured  Items  (Plate  Thicknesses,  Shape 
Dimensions,  etc.) 
o  Poisson's  Ratio 

Eleven  major  steel  producers  in  the  U.S.  were  contacted. 
The  responses  received  [34]  did  not  include  any  new  data,  but  instead  made 
reference  to  [35]  through  [39],  Generally  speaking, the  steel  producers 
indicated  that  they  do  not  collect  the  type  of  data  requested.  Manufacturing 
of  steels  is  controlled  within  the  limitations  of  References  [37]  and  [39]. 

The  literature  survey  disclosed  a  number  of  pertinent 
references  that  are  of  value.  These  references  contain  information  on 
strengths,  ^atigue,  corrosion,  dimensional  accuracy,  and  welding  stresses. 

Since  most  ship-related  probabilistic  structural  analyses 
make  use  of  only  a  few  uncertainties  and  usually  assume  values  for  their 
COV's,  some  of  the  references  cited  were  used  to  obtain  more  accurate  estimates 
of  several  COV's.  Although  there  is  not  an  abundance  of  data,  it  may  be 
possible  to  uncover  more  numerical  values  and  do  more  comprehensive  an 
analysis  than  what  is  reported  here.  This  work  should, therefore^ be  continued. 

7.2.2  Forms  of  Existing  Data 

The  data  for  the  uncer ta i nt i es  considered  in  this 
study  fell  into  three  categories: 

o  Means  and  Variances 
o  Data  points 
o  Tolerances 

For  each  of  these  a  different  approach  was  used  to  compute  the  respective  COV. 

In  the  first  case  the  COV  can  simply  be  computed  by 
dividing  the  square  root  of  the  variance  by  the  mean. 

In  the  second  case,  the  data  points  can  be  used  to 
directly  compute  the  mean  and  variance  of  the  uncertainty,  which  of  course 
would  directly  yield  the  COV: 
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p  = 


I 


N 
tJ 

£=  2 


Z  *X. 

(46) 

’l:.l 

(V.  -0* 

1  N-l 

(47) 

£ 

(48) 

where : 


N  =  number  of  data  points 
xn=  nit  data  point 
u  =  mean  of  variable  x 

o 

°x  var'ance  variable  x 


In  the  third  case,  if  it  is  assumed  that  the  uncertainties 
arise  because  of  many  individual  variables,  it  can  be  shown  by  the  "Central 
Limit  Theorem"  that  the  uncertainties  will  be  normally  distributed.  It  has 
been  determined  for  some  ship-related  uncertainties  that  the  tolerance  limits 
generally  encompass  99-7 %  of  the  events  £53] .  For  a  normally  distributed 
random  variable, this  corresponds  to  the  following: 


Pr 


P*  to] 
If  (x)dx' 
p-tol 


\  o  /2tt 
p-  to! 


(X-P) 

e*p[  20 Z  J  dx 

=0.997 


1*9) 


where :  p  (x)  =  probability  density  function  and  tol  =  tolerance  limit. 


By  the  change  of  variable  t  = 


equation  (49 )  reduces  to: 


r  o- 


Pr  =  /  7TS  «4(-t2/2)dt  =  0.997 


Cr>0) 


-isl 
cr 


Equation  (50)  represents  a  zero  mean,  normally  distributed  process  with  a 
standard  deviation  of  1.0.  It  is  known  that  the  .997  probability  level 
is  contained  within  -  3  standard  deviation.  Therefore: 


^  tol  3 

(51) 

a  = 

(52) 

It  is  mentioned  in  [54]  that  tolerance  limits  are  usually  taken  to  be 
pt3c.  Consequently,  on  the  basis  of  the  foregoing,  the  mean,  p,  can  be 
chosen  from  the  context  and  the  standard  deviation  computed  by  dividing 
the  tolerance  by  3- 


The  uncertainties  which  were  analyzed  are  listed  below 
in  groups  corresponding  to  the  type  of  data  that  were  available: 
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o  Means  and  Standard  Deviations 

-  Depth  of  stiffener  web  [53] 

-  Breadth  of  stiffener  flange  [53] 

-  Breadth  and  length  of  plate  [53] 

-  Yield  strength  [35] 

-  Tensile  strength  [35] 

-  Initial  deflections  in  plates  [67] 

o  Data  Points 

-  Depth  of  ship  [40] 

-  Beam  of  ship  [4o] 

-  Flange  breadth  [40] 

o  Tolerances 

-  Depth  of  ship  [40] 

-  Beam  of  ship  [40] 

-  Thickness  of  plate  (receipt  inspection)  [40] 

-  Thickness  of  plate  (undercut)  [40] 

7.2.3  Determination  of  COVs 


The  means  and  standard  deviations  from  [35]  are  presented 
in  Tables  I  through  3  along  with  the  computed  C0\ Is. 

In  the  case  of  the  means  and  variances  available  for 
flange  breadth,  web  depth,  arid  length  of  plate,  it  is  assumed  that  the  variance 
is  not  a  function  of  the  absolute  size  of  the  member  since  it  is  not  presented 
in  this  mariner.  Consequently,  mean  and  the  C0V  of  a  dimension  "L"  are 
given  by; 


>1  =  x+l 

(53) 

cov  =  - 

V 

(5*0 

Figures  3  through  5  give  the  data  base  and  express :ons 
for  the  uncertainties. 


The  determinations  of  CO Vs  from  data  points  for  flange 
breadth,  ship  depth,  and  ship  beam  are  shown  in  Tables  *t  through  6.  Note 
that  the  results  of  Table  and  Figure  **  agree  well  in  that  they  both  show 
C0V  around  1%  for  flange  breadth.  It  should  be  noted  that  the  COVs  for 
ship  beam  and  depth  are  extremely  small.  Reference  [10]  assumed  a  value 
of  zero  which  appears  reasonable. 

The  determination  of  COVs  from  tolerances  for  ship 
depth,  ship  beam,  thickness  of  plate  (receipt  inspection),  and  the  thickness 
of  plate  (undercut)  are  given  in  Tables  7  through  10. 

7 . 2 . COVa  From  Literature  Survey 


Reference  [11]  presents  some  results  for  objective 
uncertainties  for  materials,  scantlings,  and  manufacturing  imperfections. 
These  are  shown  in  Table  11, 
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FIGURE  5 


DETERMINATION  OF  COEFFICIENT  OF  VARIATION 
FOR  THE  UNCERTAINTIES  BREADTH  AND 
LENGTH  OF  PLATE 
(Data  from  reference  53) 
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TABLE  4 


FLANGE  BREADTH  UNCERTAINTY 


Da  ta 
Point 

x  (inch) 

(x  -  x)2 

l 

1/8 

0.000961 

2 

1/8 

0.000961 

3 

1/4 

O.OO8836 

4 

1/8 

0.000961 

z 

5/8 

0.01172 

wc=  flange  breadth 


x  *>  ~  -  0.156" 


■4 


o  =>  0.0625" 


if  N < 2 5 


Wr,  “  X  +  W(  =  0.156  +  W, 

flange  f  f 

breadth 


0.0625 


S  _ 

p  0. 156  +  w. 


For  w^.  =  6"  : 


o~:?§sVgs  0-o,og- 

Comparison  with  Figure  4 : 

Compare  2.  “  1.02%  with  that  using  formulation  from  Japanese  standard  data 


*•  6"  =»  152.4  mm  (1  inch  =»  25-40  mm) 


2.18 


»•  0.51  +  t 


flange 


irTBEr"’01*3  °r  '•« 
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TA8LE  V  -  OE TERHINATfOft  OF  6D  (COV  OF  DEPTH)  FROM  THE  MEAN  AND  STD. 

- mrw  raita  - ■ — ~ 


DEPTH  UNCERTAINTY 


Oata 

Depth 

Point 

Depth 

Hees.  Dev. 

i 

36' 

)  .5" 

2 

46* 

It  -  0.552' 

3 

16' 
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Table  8 


_ 

UNCERTAINTY  - 

BEAM  OF  SHIP 

Tol . 

( 1  rich) 

6 

( Feet) 

u  (Mean) 
(Feet) 

Coef.  of 
Variation 

I* 

<5 

(Percent) 

l 

.]% 

0.02k 

72" 

0.000333 

0.0333 

2 

1/2 

0.0139 

200 

0.0000695 

0.00695 

3 

1/2 

0.0133 

75 

0.000177 

0.0177 

4 

1/2 

0.0139 

96 

0.00011*5 

0.011*5 

6  avg 

-  .0181 

Table  9 


UNCERTAINTY  - 

THICKNESS 

(RECEIPT  INSPECTION) 

Data 

Point 

Tol. 

(Inch) 

6 

{ Inch) 

P (Mean) 

( 1 nch) 

6 

■■ 

■ 

1/8 

0.0417 

t 

0.0417/t 

4.  1  7/t 

■ 

.It 

0.0333t 

t 

0.0333 

3.33 

■ 

1/32 

0.0104 

t 

0.0104/t 

1.04/t 

■ 

1/64 

0.0052 

t 

0.0052/t 

0.52/t 

■ 

1/8 

0.0417 

t 

0.0417/ t 

4.17/t 

-41 


TABLE  10 


UNCERTAINTY  -  THICKNESS 


(UNDERCUT) 


\ 

!  Tol . 

1  (inch) 

6 

( 1 nch) 

y (Mean) 

( 1 nch) 

6 

1/32 

0.0104 

t 

0.104/t 

1/lb 

0.0625 

t 

0.0625/t 

1/32 

0.0104 

t 

0.0104/ t 

1/32 

0.104 

t 

0.0104/t 

1/32 

0.0104 

t 

0.0104/ t 

1/32 

0.0104 

t 

0.0104/t 

1/16 

0.0625 

t 

0.0625/t 

1/32 

0.0104 

t 

0.0104/t 

1/16 

0.0625 

t 

1 

0.0625/t  j 

1.0  4/t 
6.25 ft 


p 


TABLE  II:  Objective  Uncertainties 


wBimam 

Variable 

Source 

COV 

Material 

Modulus  of  Elasticity,  E 

Unknown 

2.5% 

Yield  Strengths  (Royal 

Navy  B  Steel) 

hi stograms 

(two) 

6-8% 

Scant  1 i ngs 

Plate  thickness 

(0.25"  plate) 

[55) 

3-6% 

Plate  thickness 

(2"  plate) 

[55] 

0.7% 

Plate  thickness 

(All) 

[56] 

4% 

Manufacture 

Imperfec¬ 

tions 

Residual  Welding  Stress 

Unknown 

10-15% 

Reference  [10]  gives  identical  information  on  objective 
uncertainties  as  shown  in  Table  12. 

TABLE  12:  Objective  Uncertainties 


WMnrnmm 

Variable 

Source 

COV 

Material 

.. 

Yield  Strength  (23  ksi  steel) 

157] 

6-8% 

Yield  Strength  (32.6  ksi 
mean  yield) 

[4lj 

6.7% 

Yield  Strength 

[58] 

7.9% 
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Reference  [67]  yives  information  on  initial  lateral  deflections  in 
plating  in  the  bottom  of  a  universal  bulk  carrier  and  a  tanker  as  shown  in 
Table  12-A.  It  was  shown  by  analysis  of  the  histograms  by  the  -  criteria 
that  the  initial  delfections  obey  a  Gaussian  law. 


TABLE  12-A 


RMj 

m 

a 

b 

Type  of 
Structure 

t 

(mm) 

1 

b/ 1 

B 

? 

or 

r, 

T 

UBC 

25,000 

2.0 

bottom 

18 

If6. 0 

1 .  7^0 

0.194 

0.0237 

0.15 

Tanker 

100,000 

deck 

22 

40.9 

i.5^7 

0.145 

0.0624 

Ik  43 

bottom 

24 

37.5 

1 .418 

0.136 

0.0618 

0.45 

TABLE  13 

SUBJECTIVE  UNCERTAINTIES  til} 


Mode  of  Failure 

Ship  Bending* 

Gross  Panel** 

Tension  yield . 

Plate  Buckl ing 

Strut  Panel  "1 

3 

0 

Beam  Column) 

4 

6 

Gr 5 1 lage  Buckl inq . 

4 

7.5 

*  c.o.v.  of  subjective  uncertainties  in  strength  arising  from  ship  bending 
actions.  {%) 

ft*  c.o.v.  of  subjective  uncertainties  in  strength  arising  from  gross  panel 
actions.  (%) 
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7-3  Subjective  Uncertainties 


The  subjective  uncer ta i nt ies, quot i ng  from  [3],  "include 
the  nonineasurabl e  inaccuracies  of  engineering  analyses,  any  non-measured 
variances  of  construction  and  fabrication,  and  the  unavoidable  and  non- 
de terminable  errors  associated  with  the  prediction  of  future  conditions." 

In  Reference  [10],  which  contains  ship  longitudinal 
strength  analyses,  for  tensile  yield  ductile  bending  and  inelastic 
compression  buckling,  the  COV  of  the  subjective  uncertainties  was  taken 
as  3%  and  5%  respectively.  In  the  case  of  bending,  the  uncertainty  was 
attributed  to: 

"1.  The  use  of  the  simple  beam  theory  which  is  based  on  Navier 
hypothesis.  This  hypothesis  excludes  any  shear  lag  or 
shear  deformation  effects. 

2.  The  presence  of  small  cutouts  and  openings  in  the  deck. 

3.  The  residual  stresses  due  to  welding. 

4.  The  cracks,  voids,  and  other  flaws  in  the  material." 

For  buckling,  in  addition  to  the  above,  the  uncertainty 
in  initial  deflections  of  the  plate  was  added.  No  mention  of  how  the 
uncertainties  were  specifically  determined  is  made. 

In  Reference  [ll],  subjective  uncertainties  are  given  for 
various  strengths,  including  ultimate  tension,  bending,  and  compression 
failure  due  to  plate  and  grillage  buckling.  The  discussion  therein 
cites  a  number  of  references  from  which  data  was  obtained  and  the  CCV's 
estimated.  The  results  are  presented  in  Table  13. 

7. 4  Concl us i ons 


Future  efforts,  in  the  long  and  short  terms,  should 
be  directed  towards  identifying  and  quantifying  more  uncertainties.  Most 
subjective  uncertainties  are  really  "as-yet-unquanti f ied"  objective 
uncertainties.  They  possess  the  property  of  being  measurable  although 
data  may  not  exist  or  may  be  sparse  and  difficult  to  locate. 

The  data  presented  herein  verifies  the  assumption  in  the 
literature  that  the  variability  in  principal  ship  dimensions  is  small. 

The  additional  results  for  material  yield  strength  also  agree  well  with 
what  has  been  presented  before. 
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SECT lOti  8.0 


SAMPLE  CALCULATIONS 


8.1  General 


It  was  stated  in  Section  6.3*2  that  the  "Approximate  Probabilistic 
Structural  Analysis  Method"  of  Reference  [13]  can  be  used  to  obtain  more  in¬ 
sight  Into  ship  hull  longitudinal  strength  through  determination  of  the  safety 
index,  Y • 


In  this  section,  the  method  is  applied  to  the  ten  different 
ships  in  vertical  hull  bending  mode  of  failure  which  were  analyzed  using  the 
method  found  in  Reference  [31].  Table  14,  reproduced  from  this  reference,  gives 
the  input  data  and  the  safety  index  results. 

This  section  presents  a  computer  algorithm  which  embodies  the  subject 
method.  An  equation  for  the  strength  COV  for  ductile  yield  vertical  bending 
is  developed  in  terms  of  the  COV's  of  the  uncertainties.  Results  of  the  com¬ 
putations  with  the  computer  algorithm  are  presented  for  nine  ships  using  the 
data  from  Table  14.  In  the  cas^  of  Ship  #10,  the  formula  developed  herein 
for  the  strength  COV  is  utilized  along  with  the  remaining  data  of  Table  14  to 
obtain  the  safety  index. 

8 .2  Computer  Algorithm 

The  computer  algorithm  represents  the  computational  sequence  of  Equa¬ 
tion?  (42)  thru  (44)  of  Section  6.3.2.  Within  this  f ramework, al 1  quantities  in 
the  equations  are  input.  However,  in  the  case  of  the  strength  COV, the  option 
exists  of  evaluating  it  from  the  strength  COV  equation  and  the  input  COV's  of 
the  strength  uncertainties. 

The  alqorithrn  is  structured  as  shown  in  Figure  6.  The  main  program 
consists  of  Equations  (42)  thru  (44)  noted  above.  The  subroutines  for  the 
strength  COV's  are  based  on  individual  derivation  of  the  strength  COV  determined 
by  the  technique  given  by  Equation  (26)  of  Section  3 - 3 3  and  the  strength  equa¬ 
tion.  Only  the  ductile  yielding  of  the  hull  girder  due  to  vertical  bending  is 
considered  here.  Equation  (20)  of  Section  3.3-3  is  the  strength  equation  for 
this  mode  of  failure  of  the  hull  girder.  If  the  areas  of  the  flange  and  the 
web,  Af  and  Aw,  are  respectively  expressed  as  functions  of  additional  variables, 
the  following  equations  would  result: 


Af  -  2[Btd+Mfts(ltW)  +  tfl)] 

(55) 

Aw  •=  4[Dtw+MwUW2+  2-f 2)  ] 

(56) 

S  =  [DBtd+DMftsUw,  +  Slfl)  +  (2/3)D2tw 

(571 

+  {2/3)DMwtsUw2+tf2)}sy 
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SHIP 


LBP 

(ft) 

m 

(f t-ton) 

(f t-ton) 

V 

2 

V 

s 

Y 

594.00 

206,000 

34,500 

0.1595 

0.11 

Ml 

656,20 

328,500 

42,500 

0.1400 

0.11 

5.840 

700.65 

396,500 

51,500 

0.1402 

0.11 

5.796 

719.10 

357 ,500 

48,850 

0.1444 

0.11 

5.769 

754.70 

519,500 

60,500 

0.1325 

0.11 

5.993 

775.00 

610,000 

61,000 

0.1229 

0.11 

5.856 

800.00 

613,500 

65,500 

0.1268 

0.11 

6.103 

1,000.00 

1,474,500 

113,250 

0.1105 

0.11 

6.152 

1,069.25 

1,718,300 

118,400 

0.1068  ’ 

0.11 

6.279 

1,076.00 

1,906,500 

131,400 

0.1068 

0.11 

6.314 

TABLE  14:  "APPROXIMATE  PROBABILISTIC  METHOD 
(Input  and  Results  from  Reference  [31]) 


m  =  Mean  of  Still  Water  Bending  Moment 
X  =  Average  Value  of  Long  Term  Wave  Bending  Moment 
Assumed:  mg  =  m+X 
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where  all  terms  are  as  previously  noted  and  in  addition: 


B 

D 


td  » 
ts  ■ 


*w| 

*w2 

**1 

s 

Mf,Mw 

s„ 


Beam  of  Ship 
Depth  of  Ship 
thickness  of  dec*.  plating 
thickness  of  stiffeners 

thickness  of  webs  (side  plating  on  longitudinal  bulkheads) 

length  of  deck  or  bottom  stiffener  web 

length  of  side  or  longitudinal  bulkhead  stiffener  web 

length  of  deck  and  bottom  stiffener  flange 

length  of  side  or  longitudinal  bulkhead  stiffener  flange 

number  of  stiffeners  along  deck  and  side  plating,  respectively 

tensile  yield  strength 


The  coefficient  of  variation  of  the  bending  strength  as  given  by  Equation  (57), 
is  shown  in  Appendix  A  to  be: 


«s2  ■  Inwik+m  2 1 1/9  (3A+3I»+*C«D) 

+A2(6B2+6td2)+E2(S2wi 

+F25ifl+1/9(3E+2G)24S] 


(58) 


r2  Mc2«tw+G2«L2+H25if2) 

+  6sy+  l3A+3B+2C+2D)2 


A  =  Btd 

B  -  Mfts{£w1+  £fl) 

C  “  Dt 

w 

D  "  Vs(*w2  +  t;2) 

E  0  MftsAwl 
F  -  Mftsif, 

6  “  «wts*w2 

H  "  Mwtsif2 
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No  attempt  has  been  made  to  include  strength  equations  for  hull  ultimate,  com¬ 
pression,  fatigue  and  other  time-dependent  strengths  due  to  the  limited  nature 
of  this  project.  As  mentioned  in  Section  6.C,  it  is  considered  plausible  that 
strength  equation  can  be  adopted  for  each  of  these  strengths  which  would  conform 
to  an  analysis  by  the  "Approximate  Probabilistic  Method". 

As  the  additional  strength  equations  and  the  expressions  for  strength 
CGV  are  developed  for  different  modes  of  failure,  they  can  be  added  to  the  al¬ 
gorithm  presented  herein  as  subroutines.  It  is  noted,  however,  that  the  algo¬ 
rithm  is  simple  and  could  be  coded  quickly  for  any  type  of  computing  machine. 

A  listing  of  the  algorithm  written  in  FORTRAN-IV  for  an  IBM  computer 
is  given  in  Appendix  B  along  with  the  documentation. 

S . 3  Analysis  and  Resu ]_ts 

The  computer  algorithm  was  verified  by  re-computirg  the  results  of 
Table  l4  with  the  input  data  given  therein,  including  V^,  the  strength  GOV. 

In  the  case  of  Ship  #10,  scantling  plans  were  available  and  V<-  was  alto  computed  by 
ths  subroutine  representing  Equation  (53). 

Ship  #10,  the  "UNIVERSE  IRELAND"  is  a  large  oil  tanker  which  has 
been  used  as  a  subject  in  many  research  studies.  Table  15  gives  the  prin¬ 
cipal  characteristics  of  the  "UNIVERSE  IRELAND". 

Figure  7  represents  the  approximate  midship  section  used  to  simplify 
the  application  of  equation  ([>8).  Since  the  vessel  is  full  at  midships,  the 
error  Introduced  by  this  assumption  should  not  be  significant. 

Table  i6  gives  structural  variables  determined  from  the  drawings. 


The  next  few  steps  of  the  calculation  were  actually  performed  by 
the  computer.  The  process,  however,  consists  simply  of  inserting  the 
above  values  into  Equation  (58)  to  yield  the  following: 


2  =  i.23  5p  +  0.798  +  0.7985^  +  0.189$* 


+  0.036, 


n  *  2  4-  6s  2  +  0 . 0 1  2  t  2  t  +  0.005  , 

+  0.33  ot  y  w  I 


*■0. C00l*o  lf  (59) 

'  *> 

Table  17  gives  the  uncertainty  COV's  determined  from  the  data 
presented  in  section  7.0. 


Equation  (59)  yields  $s  =  .081*9.  Assuming  subjective  uncertainties  to 
t>e  33>,  t'ne  total  strength  C0V  becomes 


ri  2 

C0V=  /6  +  A 


4^0849)  2  +  (.03)" 


=  0.09  (b0) 


It  is  interesting  to  note  the  relative  importance  of  uncertainties 
as  given  by  equations  £9  )  and  (60 ) .  First,  the  assumed  3V  C0V  for 
subjective  uncer ta i nt ies  does  not  significantly  affect  the  uncertainty 
that  would  be  obtained  from  objective  considerations  •..nly.  Table  18 
gives  the  values  of  the  individual  terms  of  equation  (59). 
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As  can  be  seen  in  the  case  of  the  ''UNIVERSE  IRELAND",  the  yield 
strength  of  the  material  far  outweighs  the  consideration  of  other  var¬ 
iables.  Additional  variables  of  some  importance  are  thickness  of  deck 
pitting  and  stiffeners  and  the  length  of  flanges  on  side  and  lonoi- 
tudinal  Dulkhead  stiffeners. 

The  results  of  the  safety  index  computation  employing  the  computer 
program  are  plotted  in  figure  8. 


TABLE  1 5 :  UNIVERSE  IRELAND  CHARACTERISTICS 


TYPE 

TANKER 

Approx i 'Male  dwt,  tons 

326,585 

0 ve ra 1 1  length,  ft- 

1,135-17 

i  ft 

BP  * 

1,076 

Breadth,  ft. 

17^.87 

Dept.  ft. 

105 

Design  draft  (keel),  ft. -in. 

~  1  #l 

81-5 

Builder 

1 shikawaj ima  Harima 

Block  coefficient  (L^) 

0.86 

Section  modulus,  top,  in^-ft. 

566,79;i 

Maximum  Stillwater  bending  moment  in  full 

load  condition  (long  voyage)  (sagging) 

=  1 ,940,000  LT-FT 

j 

Maximum  Stillwater  bending  moment  ?n  full  j 

1 

load  condition  (short  voyage;  (sagging) 

-  2,355,000  LT-FT 

Maximum  Stillwater  bending  moment  in  normal 

ballast  (hogging) 

=1,858,000  LT-FT 
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FIG.  8  -  RELIABILITY  ANALYSIS  -  DISTRIBUTION-FREE  METHOD 

(COMPARISON  OF  RESULTS  FROM  REF.  (31)  AND  FROM  COMPUTER  PROGRAM  OF  APPENDIX  B) 


SECTION  9.0 


CONCLUSIONS 


The  conclusions  reached  in  the  course  of  this  study  are  summarized 
below  under  the  major  subheadings  from  which  they  were  derived. 

General  Probabilistic  Structural  Analysis: 


o  Mechanical,  Civil,  and  Ship  probabilistic  analyses  of  structure 
have  all  utilized  similar  methods  which  differ  in  the  assumptions 
concerning  load  and  strength  distributions. 

o  Strengths  can  be  described  from  known  d i str : but  ions ,  synthesized 
from  known  distributions  of  constituent  parts,  synthesized  from 
assumed  distribution  of  constituent  parts,  or  synthesized  in  terms 
of  means  and  variances  from  those  of  the  constituent  parts.  The 
accuracy  decreases  from  the  former  to  the  latter  approaches. 

o  Since  probability  of  failure,  rargin  of  safety,  and  reliability 
involve  integration  under  the  tails  of  distribution  curves,  the 
shapes  of  the  distribution  may  be  important. 

o  Not  much  is  available  in  the  literature  concerning  time 
dependent  strength  analyses  in  general. 

Modes  of  Fa i 1 u re : 

o  The  most  significant  mode  of  failure  of  ship  hulls  from  a  longitudinal 
Strength  standpoint  is  not  known .  The  probable  modes  include  ductile 
yield  bending  of  the  hull  as  a  free-free  beam,  ultimate  plastic  collapse, 
buckling,  fatigue  and  fracture. 

o  Modes  of  failure  may  vary  according  to  ship  type,  size,  etc. 

o  Local  stresses  must  be  superimposed  on  primary  stresses  when 
considering  modes  of  failure  i.e.  principal  stresses  must  be 
known . 

o  The  lack  of  knowledge  of  ship  modes  of  failure  for  longitudinal 
Strength  point  to  the  need  to  obtain  more  insight  into  these  types 
of  failures  to  be  better  able  to  establish  guidelines  for  new  and 
structurally  different  ships  of  the  future. 

Load i ngs 

o  The  total  load  scenario  for  a  ship  at  sea  is  not  clearly  known  in 
terms  of  statistical  distributions. 

Structure  Statistics: 


o  Strength  and  uncertainty  distributions  for  ships  are  not  available 
except  for  limited  cases  of  the  latter.  Assumptions  are  made  in 
analyzing  individual  structures. 


o  Strength  statistics  in  terms  of  means  and  variations  have  been 
generated  for  ships  using  an  approximate  method  to  obtain  the 
mean  and  variance  from  those  of  the  uncertainties  in  the  strength. 

o  Little  has  been  done  with  strengths  other  than  static  bending 
strength  and  plate  buckling  strength. 

S trength  Uncertainties: 

o  It  is  difficult  to  obtain  data  on  strength  uncertainties. 

o  Dimensional  uncer ta i nt ies  in  principal  cha racter  i  s t i cs  of  ships 
appear  neg 1 i g i ble . 

O  Strength  variations  in  yield  and  tensile  properties  of  steels 
appear  significant. 

o  It  is  difficult  to  estimate  subjective  uncertainties;  not  much 
data  are  available. 

o  Uncertainty  statistics  can  be  obtained  from  tolerances  in  production, 

Shi p  Ana  1 yses : 

o  Most  ship  probabilistic  structural  analyses  have  been  concerned 
with  vertical  bending  an  i  yielding. 

o  Due  to  the  lack  of  statistics  concerning  ship  stregths  and  loads, 
the  development  of  an  analysis  approach  which  does  not  require  a 
distribution  shape  appears  warranted.  This  would,  therefore, be  a 
semi -probabi 1 i st ic  approach. 

o  Analyses  using  an  approach  of  the  latter  type  could  be  used  to 

compare  past  and  present  designs,  modes  of  failure,  etc,  ,  to  obtain 
more  insight  into  ship  longitudinal  strength. 


SECTION  10.0 


RECOMMENDATIONS 


An  overall  recommendation  is  that  investigations  should  continue 
in  the  area  of  ship  longitudinal  strength  so  that  any  structurally 
different  ship  of  the  future  can  be  properly  designed  with  confidence. 
Although  the  probabi 1 istic  structural  analysis  of  ships  will  probably 
have  to  remain  at  a  more  or  less  simplistic  or  semi -probabi 1 i st i c  level 
for  the  near  future  due  to  lack  of  methods  and  data,  the  efforts  to  im¬ 
prove  the  approach  should  nevertheless  be  undertaken  if  technology  is 
to  forge  ahead. 

Some  specific  recommendat i ons  are  presented  below  for  both  the 
short- and  long-term  goals: 

Long  Term 

°  Continue  to  develop  techniques  and  obtain  data  for 
both  load  and  strength  for  probabilistic  analysis. 

°  Perform  classical  probability  of  failure  analyses 
for  different  ships;  compare  and  update  the  results 
as  better  data  and  techniques  are  developed 

Strength  and  Load  Distributions: 

°  Determine  conclusively  whether  it  would  be  practical 
or  not  to  determine  strength  distributions  from  small 
scale  structural  models. 

°  Obtain  accurate  coefficient  of  variation  estimates 
for  strength  variables. 

°  Synthesize  strength  distributions  using  coefficients 
of  variation  of  strength  variables  and  assumed  dis¬ 
tributions  of  these  variables. 

°  Determine  the  exact  distributions  of  strength  vari¬ 
ables  and  synthesize  strength  distributions  by  the 
methods,  discussed  earlier,  to  obtain  a  more  accurate 
distribution  than  otherwise  available  with  the  assumed 
distribution  approach. 

°  Determine  whether  or  not  specific  tabulated  distribu¬ 
tions  are  accurate  for  ship  strengths. 

Strength  Equations: 

°  Develop  accurate  strength  formulas  for  hull  failure  in 
ductile  beam  bending,  compression  buckling,  ultimate 
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strength  failure,  fracture  and  fatigue. 


0  Incorporate  into  these  formulas  as  many  uncertainty 
variables  as  possible. 

Time -Dependent  Analyses: 

°  Develop  an  accurate  probability  of  failure  procedure  for 
the  analysis  of  time-variant  ship  strengths  for  the  cases 
of  fatigue,  corrosion,  etc. 

Short  Term 

°  Using  a  probabilistic  structural  analysis  method  of 
complexity  consistent  with  required  assumptions, 
analyze  past  and  current  ships,  including  those  that 
have  failed  in  longitudinal  strength,  for  various  modes 
of  failure.  This  method  should  embody  a  semi -probabi 1 i st i c 
approach.  Determine  the  corresponding  safety  factors. 

°  Compare  results  of  the  analyses  to  gain  more  insight 
into  ship  longitudinal  strength. 
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APPENDIX  A 


DERIVATION  OF  EXPRESSION  FOR  STRENGTH  COV  FOR  DUCTILE  YIELD  OF  THE  HULL  DUE  TO 
“  VERTICAL  BENDING 


Yield  Strength^  y  (A^  +  1/3  Aw)  (J) 

D=  section  depth  (between  neutral  .ixes  of  combined  plate  and 
beams  of  deck  and  bottom) 

B=  section  beam 

A^.=  total  area  of  flanges  (deck  and  bottom) 

Aw”>  total  area  of  webs  (sides  and  longitudal  bulkheads) 

S  *  yield  stress 

Take  into  account  the  following  uncertainties  in  the  above  strength  formula: 
o  06  B;  as  described  above 
o  tj  »  thickness  of  deck  plating 
o  «  thickness  of  stiffener 

o  t^f  thickness  of  webs  (side  plating  or  longitudal  bulkhead) 
o  l^j  length  of  stiffener  web 
o  1 ^ ;  length  of  stiffener  flange 

M,t  M  =  number  of  stiffeners  along  deck  and  number  of  stiffeners 
f  w 

along  side  plating  respectively. 

Af-  [Bt^  +  Mf  vlw  +  If,  )]*  2  (2) 

Aw*  [Dtw  +  Mw  (lw2  +  lf2)]'  4  (3) 


INSERTING  INTO  (l): 

S~{((  Btd  +  M#ts(lw,  +  «f|)]-2+  V3[Dtw+Mwts  (lw2+  lf2)].4) 

S*  [DBt^  +  0Mft4  (lWj+lf))  +  2/3  02tw  +  2/3  D*  tg  (lw2  +  1  f2>]  ^  (4) 


FROM  SECTION  3-3-3: 


2 

«s 


tU 


4l>2‘ 


(5) 
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isequently,  the  following  partial  derivatives  must  he  determined: 


3s  3s  3  s  as  as  as  Ss  as  ?s 


’  '28  ’  atd  ’  2<w’  at.’TT^  3if)’  5V  ’  ' 

25.  „  [Btd  +  Mfts  (Iw,  +  lf)>  V3  0tw  +  2/3  Mwts  (lw2  +  lf  )] 

D  2 

H  'DtdS  (7> 

3?  =  DB  ^  (S) 

at“d 

is  -  2/3  l)2  s.  (9) 
it 

w 


-  [0Hf  (  )Wi  +  V(  >  +  2/3  DHw  +  ^2>  ]  to 


57  -DVs^  (M) 


TTf  =0VS^  (,i) 


i_s  -  DBtd  +  DMfts(lw,  +  If,)  +  2/3  Dt  +  2/3  DM  tc(lw,  +  ]r  ) 


V,V,W2  +  V 


2/3  DM  t 


I  t  s. 

WSJ 


3T  “  2/3  Dtt-ts  ^  05);  Equation  (5)  then  becomes: 

f2 

^  "  (H  ’  f  +  (3T  *  T  >«2  +  (§f  •  |V«2td  +  (|f  Jg)2«2 

B  dtd  S  <>£w  S 

V.„ 

.  /*S  Tw  9  2  .  /is  .  V  2  3^  ./3s  *3  v2  2 

+  'iTbi  « J^\2  5  i  +  hr  — :2  )  3,f  +  (t-  •  -r  )  65 


let :  A=  Bt 


B=  M,t  (lw.+  1-  ) 
f  s  I  f  , 


E=  M  t  1 

f  5  w, 


F=  V.l 


C=  Bt 


w 


D=  Vs  (,w2 +  ]h] 


G=Vs,w2 

H*  W#2 


2  ,A  +  B  +  V3C  +  2/3  Dv2  2  ,  r  A _ v2  2 

«  s  “  'A  +  B  +  2/3C  +  2/3  D',|SD  'a  +  B  +  2/3C  +  2/3  u  4  6  B 

, _ A _ ,2  2  , _ 2C _ ,2,2 

lA  +  B  +  2/3  C  +  2/3D*  4  °  l3A  +  3B  +  2C  +  2D* 4  w 

, _ E _ v  2  -  . _ F _ .  2  2  * 

'A  +  B  +  2/3C  +  2/3D;*  +  'A  +  B  +  2/3C  +  2/3D'*6  f. 

.  i _ 2G _ , 2  2  ,  .  ,  2H _ ,2s2lf 

v3A  +  3B  +  2C  +  2D,4<S  V  1  3A  +  3B  +  2C  +  2D;*  T2 


+  .[?■  *  2/3  7.  1  .2  , 

'A  +  B  +  -2/3C  +  2>3D;  6*  t. 


+  r 


dr) 


2  -2,2 


S  "  (3A  +  3B  +  2C  +  2D)fc 


..[1/9  (  3A  +  3B  +  *C  +  2D)262d  +  A2  (6§  +  &t  )+E‘filw, 


+  D2  Is,2  +  F2  6*f  +1/9(3E  +  2G  )26ts 
l»(C262t  +  cVl  +H262  lf) 

U>  W_  T_ 


+  62*. 


(3A  +  3B  +  2C  +  2D  }‘ 


(18) 


Equation  (18)  Is  the  final  expression  for  the  strength  COV  as  a  function  of 
the  COVs  of  the  uncertainties. 
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APPENDIX  B 


LISTING  AND  DOCUMENTATION  OF  COMPUTER  PROGRAM 


INTRODUCTION 

The  computer  algorithm  was  developed  as  described  in  Section  8.0.  The 
program  has  several  features  that  should  be  noted.  Any  number  of  ships  may 
be  analyzed  for  any  number  of  modes  of  failure.  Each  ship  and  mode  of  failure 
requires  its  own  data  set,  however.  Further,  with  respect  to  modes  of  failure, 
there  is  an  option  to  input  the  strength  COV  or  have  it  computed  by  inputting 
the  uncertainty  COV's  and  computing  the  strength  COV  by  subroutines  for  that 
purpose.  Although  the  only  such  subroutine  presently  included  in  the  program 
is  that  for  ductile  yielding  of  the  hull  girder  in  vertical  bending,  the 
program  has  been  structured  to  allow  for  additional  subroutines  to  be  added 
for  other  modes  of  failure. 

Tables  B-1  and  B-2  give  the  input  cards  and  format  for  the  program  and 
Table  B-3  a  listing  which  should  be  self  explanatory. 

Table  B-A  contains  a  listing  of  the  data  cards  for  the  analysis  discussed 
in  Section  3.0,  and  Table  B-5  presents  the  computer  output. 
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TABLE  B— 1 :  INPUT  CARDS  AND  FORMAT 


CARO 

FORMAT 

NUMBER 

OF  CARDS 

SYMBOL 

NOTES 

1 

13 

1 

NS 

Number  of  Ships 

2 

13 

1  per 
ship 

NMF 

Number  of  modes  of  failure,  each  ship 
must  have  this  card  repeated.  Input 
all  data  cards  for  each  ship  before 
repeating. 

3 

4-FI0.4 

1  per 
ship 

SO,  OS, 

ZM ,  COVL 

SD  *=  Actual  Design  MOD 

DS  «=  Average  Design  Stress 

ZM  c  Mean  of  Load 

COVL  -  COV  of  Load 

I . 

13 

1 

1  FLAG 

0-lnput 

strength 

cov 

2-compute 

strength 

cov 

Flag  to  Indicate  whether  strength 

COV  will  be  input  or  computed. 

5 

13 

1 

NTS 

This  card  to  be  filled  out  only  if 
card  !|  is  "2".  This  directs  card 
given  the  direction  of  which  mode  of 
failure  Is  being  analyzed  and  hence 
one  is  needed  for  every  mode  of 
failure.  The  next  cord  should  come 
after  card  6.  Currently  only  ver¬ 
tical  bending  is  available  in  the 
program  so  that  NTS  *  I. 

f< 

1OF7.5 

1 

see  notes 

These  cards  to  be  filled  out  only  if 
card  l|  is  "2".  The  following 
variables  in  order,  as  defined  in 
Section  8.'0,  are  required: 

7 

7F 10. ^ 

1 

[*23  quan-T 

[tlties  J 

40,  6B,  6td,  *lw-|.  4lf|  .  *ts,  6Sy, 

6tw,  6 1 W2 ,  6 1 f j  •  Depth,  Beam,  td ,  ts. 

If f  *  Iwj  •  lfj»  *w»  «w-  *w? 

8 

7FI0.lt 

1 

If  ,  Sy,  SUBCO (Subject i ve 

Uncertainties  COV) 

9 

FIO. 5 

1 

COV 

This  card  to  be  filled  out  only  if 
card  A  Is  "0".  In  that  case  the 
strength  COV  is  Input  and  is  given  by 

COV  -  /6s**  &2 

6S  ■  objective  uncertainty  COV 

As  »  subjective  uncertainty  COV 
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TABLE  B-2:  INPUT 


"»  n  n  ^  •  >  a  n  .''i  i  t\  r\  n  r\  r\  r\  n  r, 


TABLE  B-3:  LISTING  OF  COMPUTER  PROGRAM 


OMPUTER  ALGORITHM  TO  DETERMINE  SHIP  HULL  LONGITUDINAL  STRENGTH  SAFtTY 
INDICES  BY  THE  PROBABILISTIC  METHOD  DISCUSSED  IN  SR-241 

READ  IN  NUMBER  OF  SHIPS  TO  BE  CONSIDERED 


REAL  .MF.LWl  |LF1  .X.il.LW?  *LT2  .MSI 
0 1  MENS  I  ON  AUO’ 

N  »  2 
*  3 

RP  AD  (N.l)  NS 
DO  LOOP  FOR  NUMBER  OF  SHIPS 
DO  200  N3FL«1.NS 

READ  NUMBER  OF  MODES  OF  FAILURE  TO  3£  CONSIDERED  FOR  SHIP  ♦  I • 
RFADIN.1I  NV.F 

1  FORMAT  1131 

DO  LOOP  FOR  MODES  OF  FAILURES 
DO  100  I«1.NMF 

READ  ACTUAL  DESIGN  STRENGTH. SO .AVERAGE  DESIGN  FAILURE  STRESS .OS «  Ml  AN 
OF  THE  LOAD . ZM« AND  COV  OF  LOAO.COVL 

READ  (N.2 ISD.DS.2M.CCVL 

2  FORMAT  (4F10.4I 

SELECT  BETWEEN  INPUTTING  STRENGTH  COV  OR  COMPUTING  IT  BY  SUBROUTINES 
AND  INPUT  UNCERAJNTY  COVS  AND  ONE  VALUE  FOR  SUBJECTIVE  COVS 

READ  (N.l) (FLAG 
IF  ( 1- I  FLAG  I  3.3.4 

3  CONTINUE 

C  READ  IN  TYPE  OF  STRENGTH  TO  BE  CONSIDERED 
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n  o  ^  /•»  n  n  >  n  imi  n  n  n  ^  n  (■> 


TABLE  B-3:  LISTING  OF  COMPUTER  PROGRAM  (CONT.) 


°age  2 

RFAD  (N.l)  NTS 
SO  TO  ( 5  I .NTS 
*>  CONTINUE 

RC40  IN  VERTICAL  PENDING  STRENGTH  UiNCERNTAJNTY  COV ' S 

READ  (N.6)  A 
fi  FORMAT  (10F7.5I 

RF AO  IN  SHI?  PARAMETERS  NEF.9E0  TO  COMPUTE  STRENGTH  COV 

RE  AO  (,n>20)OEPTHiREAM«TO*TS>MF  »LN1*LF1  fTW»Mrf«LV<2  .LT2.SUBCC 
20  FORMAT  I7F10.4) 

CALC  SUBROUTINE  TO  COMPUTE  STRENGTH  COV  FROM  UNCERTAINTY  COV'S 

CAUL  VERTS  ( A *0E° TH. REAM. TO .T S tMF ,L'Wl ,LF 1  * TW > AM «UW2 ,LT 2 *CuVU ) 
C0v*50RT I COVO**2+SURCO**2 I 
SO  TO  7 
4  CONTINUE 

REAR  ThE  GIVEN  STRENGTH  COV 

READ  (N.fll  COV 
3  FORMAT  (F10.SI 
7  CONTINUE 


COMPUTE  "EAN  OF  STRCNoTH 


MSI *SO»OS 
THETA=msI/£M 

COMPUTE  The  safety  I OE x 

SI  =  ( THETA- 1.0  I /< SORT ( THE TA**2»COV**2 ♦COVL**2 )  I 
OUTPUT 

WRITE  (V* 10)  NSFL 

10  format  (IX. 'SHIP  NUMBER  *  'll?) 

’  I  T  E  (■■■,11)  NMF 

1  l  format  (U.'VODt  Oc  FAILURE  UU.'.MER  *’  .13) 

■/RITE  ( ■• ,  1  ?  )  S3 

17  FORMAT  (LX. 'ACTUAL  Or SIGN  MODULUS  =’>F20.4) 

WRITE  (Mi  13)  03 

13  Fo’  (AT  (IX, 'AVERAGE  DESIGN  FAILURE  STRESS  *'.F20.4I 

<M,l4t  (.’SI 

14  F  T  (IX. 'MEAN  OF  STRENGTH  ■  '.F20.4) 

V  I  if.  (M.15)  2 M 

15  FORMAT  (IX, 'MEAN  OF  LOAD  •  '  .F2U.4) 

WRITE  (M.l  ,1  COV 

1A  FORMAT  (IX, 'COV  OF  STRENGTH  «  '.F10.4) 

WRITE  IM.17)  COVL 

17  FORMAT  (IX. 'COV  OF  LOAD  *  ’.F10.4) 

•VR  I  T  E  ( .M .  1 9  I  SI 

14  FORMAT  (IX. 'SAFETY  INDEX  »  ',F1'J.4///I 

10?  CONTINUE 
20C  CONTINUE 
call  exit 


■J  o  u  u 


SU9R0UTINE  T 0  COMPUTE  THE  DUCTILE  YIELDING  STRENGTH  COV  FOR  A  SHIP  IN 
VERTICAL  9*  NO  I NG.FORmULAT I  ON  ACCORDING  TO  THAT  GIVEN  IN  SK-241. 

SU9ROU1 INE  VERTS  ( A .OEP Th «6t AM . TD i T S *MF . LW1 *LF 1 . TW *  MW i LW2 iLT i • 
1XOVI 

REAL  MF,LW1«LF1,MR.LW2iLT2 

DIMENSION  AI10I 

A 1  *9F  AH*  TO 

9«rF*TS*(LRl  +  LFl  I 

C«DEPTH*TW 

r>»-1'V«TS*  (  LW2*LT2  I 

E«,4F«TS*LWl 

F»VF*T  S*LF  I 

G«mw*  T  S«L  W2 

H»'-1W*  T  S*L  T  2 

COVl*  <  3.0»Al  +  3.0»9-*-4.U*C+2.0*Dl«*2*A(  1  )»«2/9.G 
COv2«Al**2»(A(21*»2*A(3)**21 
COV3*E**2*A(4)«*2 
COV4  *F  »  *2  » A  I  5  I  *  *  2 

COV5M  3.0«E*2.0*GI**2*A(6I**2/9.0 
COV6=9.0/ ( 3.0*A1+3.U*P+2.0*C+2»C*D>»*2 
COv7*COV5*(COVl+COV2+COV3+COV4KOV3 I 
COV3»A(  7 ) »»2 

C0V9*  4  •  0  »  (  C»»2*AI<il«»2+G»*2«A!9  )  «»2+H**2»A  (101**21 

COV 10= ( 3.0»A1*3.0»9»2  t0*C  +  2.0*D)*»2 

C0V11=C0V9/C0V10 

XOV  =  SORT  (COV7+COV8  +  COVU  ) 

RETURN 

END 
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TABLE  B-4:  LISTING  OF  DATA  CARDS 


’/  xfio 

l? 


61200. 

A 

t  •  1 

l 

lf.2. 

240600. 

-i5*? 

81323. 

3 

11 

1 

15.2 

3?1OO0. 

‘HfJOO 

•  7  300. 

o' 

.11 

1 

15.2 

44f603. 

•UtiU 

87953. 

0 

11 

1 

15.2 

435353. 

•iff* 

131533. 

C 

11 

1 

15.2 

530033. 

.1325 

143303. 

0 

11 

1 

15.2 

671333. 

.1229 

15*300. 

n 

ll 

1 

15.2 

679300. 

.1256 

355333. 

0 

11 

1 

15.2 

1597750. 

.1135 

433633. 

0 

11 

1 

15.2 

1936733. 

.1355 

433750. 

0 

11 

1 

15.2 

2337930. 

.1265 

433753. 

3 

ll 

1 

15.2 

2456430. 

.1365 

4*3750, 

7 

1 

15.2 

2<"4S  400  • 

.1253 

.331 

3015  .0332  .333 

.393  .349  .3750 

.253  .3036 

37.0 

53.3 

.035 

.225  50. 

.532 

.02 

30. 

•  5  5 

.133  .03 
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TABLE  B-5:  COMPUTER  OUTPUT 


SHIP  NUMBER  «  1 

VODF  OF  FAILURE  NUMBER  ■  1 

ACTUAL  DESIGN  MODULUS  *  51200.0079 

AV  -  RA'jF  D'SIGN  FAILURE  STRESS  *  15 

MEAN  OF  STRENGTH  •  778240.0017 

VS AN  OF  LOAO  «  240500.0316 

COV  OF  STRENGTH  •  0.1100 

CO V  OF  LOAD  *  0.1595 

SAFETY  INDEX  *  5.7123 


SHIP  NUMBER  *  2 

•<00F  OF  FAILURE  NUMBER  *  1 

ACTUAL  0ES1GN  MODULUS  «  81320.0156 

AVERAGE  DESIGN  failure  STRESS  «  15 

VFA'I  OF  STRENGTH  ■  1236J64.0034 

VEAN  OF  LOAD  *  371000,0632 

COV  OF  STRENGTH  *  0.1100 

COv  OF  LOAO  *  0.1400 

SAFETY  INDEX  *  5.9434 


SHIP  NUMRFR  *  3 

“'ODE  OF  FAILURE  NUMBER  «  1 

ACTUAL  DESIGN  MOOULUS  *  97000. 015U 

AVERAGE  DESIGN  FAILURE  STRESS  *  15 

MEAN  OF  STRENGTH  *  1 474 A 00 . 00 34 

’'FAN  OF  LOAD  »  -  448000.0632 

COV  OF  STRENGTH  *  0.1100 

COV  OF  LOAC  =  0.1402 

SAFETY. INDEX  =  5.9015 


SHIP  NUMBER  *  4 

MODE  OF  FAILURE  NUMBER  *  1 

ACTUAL  DESIGN  MOOULUS  «  87950.0156 

AVERAGE  DESIGN  FAILURE  STRESS  *  15 

■‘FAN  OF  STRENGTH  *  1 335840. 0029 

•■■CAN  OF  LOAD  =  406350.0633 

COV  OF  STRENGTH  *  0.1100 

COV  OF  LOAD  a  0.1444 


-75- 


V 


Am.' 


2000 


2000 


2000 


2300 


TABLE  B-5:  COMPUTER  OUTPUT  (CONT.) 


SAFETY  INDEX  ■  5.9770 


SHIP  NUMBER  «  5 

MODE  OF  FAILURE  NUMPER  •  1 

ACTUAL  DESIGN  MODULUS  *  131600.0316 

AVERAGE  DESIGN  FAILURE  STRESS  ■  15.2000 

MEAN  OF  STRENGTH  -  2000320.0034 

MEAN  OF  LOAD  *  590000.1264 

COV  OF  STRENGTH  «  0.1100 

CO'/  OF  LOAO  *  0.1325 

SAFETY  INDEX  «  6.0939 


SHIP  NU'-'PEP  a  (, 

MODE  OF  FAILURE  NUMBER  »  1 

ACTUAL  DESIGN  MODULUS  «  143000.0316 

AVFRAGE  DESIGN  FAILURE  STRESS  ■  15.2000 

MEAN  OF  STRENGTH  •  2173600.0058 

MEAN  OF  LOAD  *  671000.1267 

COV  or  STRENGTH  *  0.1100 

COV  OF  LOAC  «  0.1229 

SAFFTY  INDEX  *  5.9410 


SH I °  MUMPER  a  7 

'•'ODE  OF  FAILURE  NUMRER  *  1 

ACTUAL  DESIGN  MODULUS  *  15H000.0316 

AVERAGE  DESIGN  FAILURE  STRESS  =  15.2C0U 

MEAN  of  STRENGTH  «  2401600.0068 

v.;AN  OF  LOAO  *  679000.1267 

COV  OF  STRENGTH  *  0.1100 

COV  OF  LOAC  *  0.1263 

SAFETY  INDEX  =  6.1997 


SHI°  NUMPER  c  3 

'••ODE  OF  FAILURE  numrer  =  1 

ACTUAL  DESIGN  MODULUS  =  366300.0633 

AVERAGE  OFSION  FAILURE  STRESS  •  15.2000 

MEAN  OF  STRENGTH  »  5567760.0117 
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TABLE  B-5:  COMPUTER  OUTPUT  (CONT.) 


■'FAN  OF  LOAD  a 

CUV  or  STRENGTH  •  O.llGO 

cov  or  load  »  o.l 10S 

Safety  INDEX  *  6.2471 


SHI  3  NUMBER  «  9 

MODE  or  FAILURE  NUMBER  »  1 

ACTUAL  DESIGN  MODULUS  *  439600. "63? 

4VERAGF  DESIGN  FAILURE  STRESS  »  13.2000 

"FAN  OF  STRENGTH  *  6631920.0117 

MEAN  OF  LOAD  *  1836700.2334 

COV  Oe  STRENGTH  »  0.1100 

COV  OF  LOAD  «  0.1063 

SAFETY  INDEX  *  6.3691 


SH 1 °  NUMBER  *  ID 

'■'OOF  OF  FAILURE  NUMBER  *  1 

ACTUAL  DESIGN  MODULUS  ■  493750.0633 

AV-RAGc  DESIGN  FAILURE  STRESS  •  15.2000 

"E  AN  OF  STRENGTH  «  7505000.0136 

VTA'!  OF  LOAD  *  2037900.2529 

COV  OF  STRENGTH  *  0.1100 

cov  OF  LOAD  *  0.1063 

$ ArcT Y  INDEX  *  6.4035 


SHIP  NUMBER  *  11 

vODE  OF  FAILURE  NUMBER  =  1 

ACTUAL  DESIGN  MODULUS  *  493750. C633 

AVERAGE  DESIGN  FAILURE  STRESS  =*  IS.  2000 

"FAN  OF  STRENGTH  *  7505000.0136 

MEAN  OF  LOAD  *  24  96400. 505a 

COV  OF  STRENGTH  »  0.1100 

Coy  OF  LOAD  *  0.1058 

SAFETY  INDEX  *  5.7370 


SHIP  NU'I’FR  >  12 

voor  OF  FAILURE  NUMBER  *  1 


AC  <UAL  DESIGN  MODULUS  »  492750.06  53 

AVERAGE  DESIGN  FAILURE  STRESS  a  15.2000 

MEAN  OF  STRENGTH  «  7505000.0136 

' * r A N  OF  LOAD  «  24*6400. 5058 

COV  OF  STRENGTH  «  0.0340 

COV  OF  LOAD  •  0.1068 

Safety  INDEX  «  7.3332 


•O.S.  OOVOSOtEHT  PBIKTIHG  OmCE:  198 1-0-725-969/' 5. 1 
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